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Efficiency of 47r—Crystal-Scintillation Counting: 
1. Experimental Technique and Results 


C. C. Smith, H. H. Seliger, and J. Steyn ' 


The technique of 4z-crystal-scintillation counting has been applied to the standardiza- 


tion of beta-emitting nuclides. 


Phototubes viewing an anthracene 


sandwiched source at 


180° have been used with high-gain nonoverloading amplifiers in addition and in coincidence 


circuits to obtain high detection 


analyzer is described. The 
Tl are 


1. Introduction 


The method of 478-gas-proportional counting has 
been established for some 6 years now for the primary 
standardization of beta- emitting nuclides at the 
National Bureau of Standards [1,2 2,3)" The method 
of 4z-liquid-scintillation counting has also recently 
been investigated at the Royal Cancer Hospital in 
London [4]. It was therefore considered to be of 
interest to initiate experiments to determine the 
efficiency of 4x-crystal-scintillation counting and its 
applic ability to primary beta st: andardization. 

These measurements have now been concluded 
and show that the method of 42-erystal-scintillation 
counting, although not as generally applicable as 
4r8-gas-proportional counting, can be used for 
medium and high-energy beta-emitting nuclides. 
It is the purpose of this paper to report briefly the 
experimental results obtained for the benefit of other 
interested workers in the field. 


2. Experimental Procedure 


The crystal scintillation counters used in these 
experiments consisted of sandwiches formed by two 
transparent anthracene crystals each approximately 
lin. in diameter and ranging from \ to % in. in 
thickness. Two Dumont 6292 multiplier photo- 
tubes, in contact with the backfaces of the anthracene 
crystals, were used with high-gain nonoverloading 
amplifiers to detect light pulses occurring within the 
crystals. One of the experimental arrangements is 
shown in figure 1. The phototubes were inside a 
light-tight refrigerator, with the lowered tempera- 
tures being used to reduce the random thermal 
tube noise. The entire sandwich is surrounded by 
an aluminum cylinder for complete light collection. 
In this connection it was found that a 2-in. diameter 
cylinder gave slightly better light collection than one 
of only slightly larger diameter than that of the 
crystals, presumably because the former enabled 
reflec ‘ted light to reach the phototube directly without 


‘Guest worker from the National Physical Laboratory, Pretoria, Union o 
South Africa, from January to June 1953 


? Figures in brackets indicate the lit rature references at the end of this paper. 
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efficiencies. 
results of measurements made with P®, Co, 
compared with those obtained by 428-gas-proportional counting. 


delay-line coincidence 
Sr” Y*, J'31, and 


A precise shorted 


being retransmitted through the crystals or scattered 
by them. 

Negative high voltage was supplied to the photo- 
tubes so that their anodes could be connected directly 
to the cathode-follower input grids. Voltages as 
high as —1,400 v were used and electronic amplifier 
gains of 40,000 were, in some cases, found necessary. 

The radioactive material was generally introduced 
into the sandwich by depositing an aliquot of the 
active solution directly onto the top face of the lower 
crystal. The deposit was allowed to dry and then 
covered by the second crystal to form the sandwich. 
Where special techniques were used to precipitate 
the active material the sources were prepared on 
thin Formvar-polystyrene laminated films (10ug/em?) 
and then placed on the crystal face. 

The pulses from the phototubes arising from the 
light pulses in the crystal sandwich were counted in 
two ways, namely in addition and in coincidence. 
In the former method the separate phototube- 
cathode-follower output pulses were added together 
and then fed into a single amplifier and counted. It 
was found that a single tube and hemispherical re- 
flector arrangement was not as efficient as the two. 
tube addition method. In the coincidence method 
the output pulses were fed into separate amplifiers, 
thence into a coincidence analyzer of short (lysec) 
resolving time and the coincidences counted. 

Each of these methods, two-tube addition and 
coincidence, has advantages and disadvantages. In 
the coincidence method the random noise pulses origi- 
nating in each phototube will not be counted except 
for those few accidental coincidences arising during 
the finite resolving time of the coincidence circuit. 
However, it is necessary that a true light pulse 
originating inside the sandwich be “seen” by both 
phototubes. This stringent requirement is absent 
in the addition method so that the sensitivity to 
very small light pulses is higher by at least a factor 
of 2 as compared with the coincidence method. On 
the other hand, the addition method suffers from 
the limitation of excessive phototube noise at the 
high tube gains sometimes required. Because the 
phototube noise is usually not reproducible to better 
than +5 to 10 percent the uncertainty in the knowl- 
edge of the true disintegration rate can become 
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appreciable. A noise counting rate of 2,000 counts 
per second in each channel would, in the addition 
method, give rise to approximately 4,000 counts per 
second to be subtracted from the total 
counting rate. However, in the coincidence method, 
with a resolving time of 1 usec, the increase in coin- 
cidence counting-rate observed would be only 4 
counts per second. Because this number is a small 
fraction of the “‘true”’ counting rate any uncertainty 
in the knowledge of the exact noise background will 
be reduced in this ratio. 

The coincidence method presupposes that at least 
one photoelectron is emitted from each photocathode 
in order that there be a finite probability of detecting 
a coincidence. Assuming a5 percent photocathode 
photoelectric efficiency and a geometrical light- 
collecting efficiency of 100 percent, one would require 
a minimum of 40 photons (20 photons to each photo- 
tube) to give rise to a coincidence count. However, 
in order to insure that both photocathodes receive at 
least 20 photons with a 99 percent probability it 
would require a minimum of 68 initial photons. 
With the further assumption of a quantum efficiency 
of 65 ev/photon for anthracene [5] the effective beta- 
ray cutoff energy would correspond to 4.4 kev. In 
the addition method, however, either phototube 
would have a threshold of 20 photons. This require- 
ment is not as stringent as in the coincidence method 
and the effective cutoff energy would be 2.2 kev. In 
the 478-gas-proportional counting method the effec- 
tive energy cutoff is estimated to be only 0.2 kev. 

For high-energy beta emitters such as P® and 
sr”-Y” the preponderance of the emitted beta rays 
have energies greater than 2.2 kev. Further, the 
light pulses are rather large so that the phototube 


observed 
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gain need not be high. Thus the noise is low and the | 
addition method is satisfactory. However, for the 
lower-energy beta emitters the phototube gain must | 
be increased, increasing the random tube noise and 
necessitating the use of the coincidence method. 

A precision shorted delay-line coincidence analyzer 
was designed for the coincidence measurements and 
is shown schematically in figure 2. The unit accepts 
positive pulses from the nonoverloading amplifiers, 
shapes them and feeds the shaped pulses into a biased 
mixing stage. The rise and decay times of th 
square-shaped pulses are 0.06 usec. The resolving 
time is variable by means of a shorted delay-line and 
the resolving time is independent of counting rate up 
to 5,000 counts per second in each channel. The 
single channel outputs record only those pulses that 
have been accepted by the input discriminator and 
which are on their way to the coincidence mixing 


stage. 





3. Results and Discussion | 


A good fraction of the data was taken at room 
temperature. Operation at temperatures below 0° C 
considerably improved the precision of the tube-noise 
measurements and permitted the use of higher 
phototube voltages than was previously feasible. 


To illustrate the improvement due to lowered 
temperatures and consequent lowered tube noise, 
as well as the difference between the coincidence 


method and the addition method, several relative 
efficiency curves have been reproduced in figure 3. 
The efficiencies relative to 1 r8-gas-proportional 
counting are plotted as functions of phototube 
voltage for both the addition and the coincidence 
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Figure 3. Relative efficiency of 4x-crystal-scintillation count- 


ing to 478-gas-proportional counting. 


methods at 28° C and at 6° C. In each case 
the usable plateau is considerably lengthened at the 
lowered temperature. These data were obtained for 

i high-energy beta emitter, namely P*. For low- 
energy beta emitters the plateau shrinks to a very 
small region of phototube voltage and the improve- 
ment due to lowered temperatures even more 
striking. For example, TP’ (2ma,=0.67 Mev) 
exhibited a plateau only at the lowered temperature 
for the addition method. At room temperature the 
thermal-noise background at the phototube gains 
required was so large and erratic as to completely 
mask any plateau. Because of its better noise- 
discrimination characteristics the coincidence method 
was used for all of the lower-energy beta emitters, 
at room temperatures and at the lowered tempera- 


is 


tures. It is particularly in these cases that the cor- 
rections for dead-time and coincidence resolving 
time become important. In a separate paper [6] 


Shorted delay-line coincidence-analyzer circuit. 


both approximate and rigorous derivations for these 
counting losses are given and their range of applica- 
bility is discussed. Because it is shown that the 
approximate derivation can be employed under the 
present experimental conditions, only these correc- 
tions will be summarized here. 

In channel 1, consisting of phototube I, amplifier 
| and discriminator and gate I, a counting rate N(1) 
is observed with 


N(I)=Nit+ Mee. (1) 


Ni is the observed thermal tube noise of phototube I 
and N,,. is the observed “true” counting rate due 


to the source. Similarly in channel II for photo- 
tube II viewing the opposite side of the sandwich 


(2) 


The dead-time losses are due to the introduction of 
an electronic gate rt, (t~1=Tgn) Whose duration is in 
all cases longer than the resolving time of the non- 
overloading amplifier. The use of this fixed gate 
insures a precise correction for dead-time as the 
amplifier resolving time is dependent upon the 
degree of overload or gain. To a first approxima- 
tion the fractional dead-time loss due to noise and 
“true” counts in both channels is given by 


(ID=Ne+ Me. 


L=7, (ND+NUID—Niet+NiNet,), (3) 
remembering that .V,. occurs in both channels simul- 
taneously and should be corrected for only once. 
Again to a first approximation the observed coinci- 
dence rate N, is related to N,. by 


N.=Nie+27-NN;, (4) 
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where +r, is the resolving time of the coincidence 
analyzer; the last term in eq (4) being the accidental 
coincidence rate due to random noise. From eq (3) 
and (4) one obtains the total disintegration rate, 


N,, as 


N.—2reNLN. 
1—7,|N(1)+NID—(N.—2teNaNw)+NLN- 7, 
(5) 


N, 


In a blank experiment, i. e., no source present, 


(N.)n= (Neen t+2re(NNo)p (6) 

CN,.)») arises from cosmic rays, possible contamina- 
tion of the erystals, and fluorescence in the glass of 
the phototubes. In actual practice (V,.))»><(NU,LD)» 
so that with very little error Ni in eq (5) and (6) can 
be set equal to (N(L) dp. Likewise N’’~(N(ID) py. 


The disintegration rate of the source N, is then 
given to within a few tenths of a percent by 
No N, ( Ne lb (7) 


Figure 4 shows the extent of these corrections in the 
cases of TP’ and Co”. The dashed curves represent 
the observed values of \V.—CN,), as functions of 
phototube voltage for both isotopes at 28° C and at 

6° C. The solid curves, representing these data 
corrected by means of eq (5), (6), and (7), give No, 
the actual source disintegration rate. Due to the 
high noise counting rates at room temperature the 
corrections required were appreciably higher than 
those for the lower temperature, making it imprac- 
ticable to continue any measurements at room tem- 
perature past 1,200 v. 

Table 1 gives a summary of data obtained by 
means of the 4x-crystal-scintillation technique for the 
nuclides P*, Co™, Sr”®— Y”, L, and TP", compared 
with measurements made in the 428-gas-propor- 
tional counter. 

This comparison points up a serious limitation of 
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solute disintegration rates for sources of cobalt-60 and thalliunm- 
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the 4x-crystal technique for low-energy beta emitters 
Because the minimum pulse observable is that caused 
by a single photoelectron being emitted from the 
photocathode there is a definite energy cutoff that 
will depend upon (a) the intrinsic photon efficieney 
of the crystal, (b) the transmission of the crystal to 
its own light, (c) the external light-collecting geom. 
etry, and (d) the photoelectric efficiency of the 
photocathode. The variations in observed efficiency 
for Co” shown in table 1 are probably due to com. 
binations of these factors. 

It is probably the case that for low-energy beta 
particles the photon efficiency of crystal anthracene 
will be lower than 1 photon per 65 ev [5]. If one 
assumes an average efficiency of 95 percent for the 
4m sandwich-crystal method in the case of Co® jt is 
possible to make some estimate of this value for 
anthracene for low-energy beta particles. Ninety- 
five-percent coincidence efficiency corresponds to 97.5- 
percent single-channel efficiency or a 2.5-percent logs 
of low-energy beta particles. As roughly 2.5 percent 
of the Co” beta spectrum is below 5 kev in energy 
this would correspond to a minimum energy cutoff 
of 5 kev. Again assuming a 5-percent photoelectric 
efficiency for the phototube photocathode one ob- 
tains a photon efficiency for anthracene of the order 
of 1 photon per 120 ev in the range 0 to 5 kev. This 
is not an unreasonable figure, for in this range the 
specific ionization of electrons is a very rapidly rising 
function in the direction of decreasing energy. 


TABLI ] Summary 4 4 f data for various nuclides 
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4. Conclusions 


The 4x-crystal-seintillation technique is suitable 
for the standardization of medium and high-energy 
beta emitters and in this energy range gives results 
that agree quite well with results obtained with the 
4r8-gas-proportional counter. One possible ad- 
vantage of the former method is the higher density 
of anthracene relative to the counter gas and _ the 
consequently higher efficiency for the detection of 
X-rays from electron-capturing nuclides. This may 
prove useful in the standardization of electron- 
capturing nuclides. 
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Efficiency of 47-Crystal-Scintillation Counting: 


2. Dead-Time and Coincidence Corrections 


W. B. Mann and H. H. Seliger 


The dead-time and 
counting have been both approximately 
under special conditions the 
additional correction terms 


requires 


1. Introduction 


using the 42-crystal-sandwich- 
scintillation counter, initiated by one of us [1]! for 
the standardization of radioactive sources, have been 
described, namely, the “addition” and ‘‘coincidence”’ 
methods. The former may be used where the counts 
due to noise in the multiplier phototubes are low 
compared with the counts due to the radioactive 
source, Whereas the latter method is applied where 
the “true”? counts due to the source are low com- 
pared with the noise counts and it is desired to 
reduce the effect of these random noise counts. 

In the addition method the dead-time correction 
is completely orthodox, it being necessary only to 
allow for the true counts lost due to the total dead 
time arising from a given number of counts in the 
one channel, each count having associated with it 
agate or dead time r¢. 


Two methods of 


2. Approximate Derivation 


In the coincidence method, however, a non-random 
element is introduced that makes the calculation of 
the dead-time and coincidence losses somewhat more 


dificult. But even so, provided that the counting- 
rate losses due to both noise and source, or true, 
counts are low, an approximate approach to the 


problem may be used. 

Let us therefore assume, in the first place, that in 
the coincidence method, the counting rates due to 
noise in the two channels (V,{ and NY per second) 
and that due to the radioactive source (NV; per second) 
are such that there is no appreciable interaction. 
If, further, Ni. is the number of true counts per 
second traversing each channel and being recorded 
as coincidences, while N(1) and N(IT) are the total 
number of counts observed in each channel respec- 
tively, then 


(1) 


Figures in brackets indicate the literature references at the end of this p 


per 


coincidence corrections 
and rigorously 
non-randomness due 


to be applied in 42-erystal-scintillation 


derived. It has been shown that 


to “‘true’’ counts appearing in both channels 


and 


N, 


NAD + Nic (2) 
If \,. were extremely small, then the dead times in 
ach channel for a gate dead time of 7, would be Nir, 
and Nats; and the probabilities of detecting a true 
count in channels I and IIT would be, respectively, 
P, —Nnatg) (3) 
and 
” 
Pr (1 ~Natg)- (4) 
As these noise counts are completely random, the 
probability of a true count being observed in both 
channel I and channel II, that is to say as a recorded 
coincidence, is given by the multiplication together 
of the probabilities given in eq (3) and (4). That is, 


Fa { ] 


a ” rv? Tw" 9 7 
1—(Ni+Na)te+NaNaz7?, (5) 
the last term on the right-hand side representing the 
overcorrection on account of overlapping of noise 
counts in the two channels. 

If, on the other hand, there were negligible noise 
counts but a number of true counts N.. per second, 
then the dead-time correction would be represented 
by a probability Pt, where 

P$=1—Neetg. (6) 

When the quantities Nj7, and Nir, are small, as in 
the conditions of the present experiments, it can be 
assumed that there is no appreciable interaction be- 
tween noise and true counts, and that the total prob- 
ability of observing a true count as a coincidence ; is 
given by the product of the probabilities of eq (3), 
(4), and (6), 

, ” r ad ww" 4 ol 

P.~l (NatNat Nee) tet NaNate, (7) 
neglecting N,.Ni72 and N,.Ni72, which are usually 
— J 

an order of magnitude less than NaNove_ NINt72 
itself involves a correction of at most 0.5 percent, 


even at high noise counting rates. 
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Substituting for N, and Nf from eq (1) and (2), 

we have 
P.~1—(N()D+ NUD —Nie) t2+- Na Nat? (8) 

and 

Nvc= Ni {1—(N(D 4+ NID — Mi) te + NaNnzz}. (9) 
Equation (8) neglects an overcorrection of the form 
(7/2) Ni+Ni), where rv, is the resolving time of the 
coincidence analyzer. However, as 7r.<r, and we 
have already essumed that Nyry and Nir, are small, 
the approximation should still be a good one 

The coincidence count will also include, however, 
accidental coincidences arising from the random 
noise counts. These accidental coincidences are 
equal to 27.N,NZ per second, where +, is the resolving 
time of the coincidence enalyzer. Therefore, the 
observed coincidence-counting rate of NV. per second 
is related to the true-count coincidence rate by the 
expression 


N.=N,.+21r.N,N,.- (10) 
From eq (9) 


Ni, 


N : ; saw, (11) 
© 1-—(N(D+N(ID—Nie) tg + Na Nar? 

N,,. is not directly observable because the observed 
coincidence-counting rate, .V,, is the sum of N,. and 
the accidental rate due to noise. Substituting for 


N,. from eq (10) into eq (11) we have that 


N.—2r.N,Nz 
1—{N(I)+N(UID—(N.—27-NiNe} e+ NiNo 32 


(12) 


N, 


In a blank experiment (without a source in the 
erystal-sandwich scintillator) 


(Neh (Neen T 27.(Na)o( Nn): (13) 


In this expression, (.N,.), will arise from cosmic 
radiation, possible contamination of the crystal 
faces, and fluorescence in the glass of the multiplier 
phototubes due to potessium-40. However, in the 
individual channels, (N,.),<.N(1),. N(ID),. so that 
within a few tenths of a percent, 


Na=(Na)o=(N(D)o— (Neln +27e(N(D NUD») 


and 


Na=(Na)e=(N(UID)p—(Nedp +27-(N(D NII 


Substitution of these values for NV; and N? in ed (12) 
and subtraction of (N,,)), as given by eq (13), from 
N,, as calculated from eq (12), vields the true dis- 
integration rate of the source to within a few tenths 
of a percent. 
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This approximate derivation of the dead-time 
losses is, however, based on certain premises that 
may not always be strictly applicable in practiee 
In the first place, if the noise counts are negligible 
the dead-time losses would be given ex xactly by ex 
(6); but in practice the noise counts can never hp 
neglected. If, on the other hand, the true counts 
are very low, the correction for dead-time logge 
would be given by eq (5); but then the method 
would be impractical, as long periods of counting 
would be required to give good statisties. The third 
possibility that has been assumed in the derivation, 
namely, that of reasonable numbers of noise and 
true counts but limited interaction between the two 
due to low counting rates, reflects the conditions of 
the present experiments. 

It was felt that several of the : assumptions made 
in the above approximate derivation are in a measure 
intuitive, and that the experimental results obtained 
would be placed on a more firm footing if a rigorous 
derivation were developed to establish the relative 
orders of magnitude of the various interactions, — [p 
addition, once one enters into the realm of very high 
counting rates with large counting-rate losses, there 
will be an interaction between noise and true counts, 
and it will be impossible to multiply the probabilities 
Pe and Pt to give the combined probability of 
observing a true count coincidence, 2., of eq 
Moreover, because a great many true counts : 
channel I are linked to true counts in channel IT, i 
will not be permissible to derive the probability * 
multiplying together two individual probabilities for 
the separate channels as was done to derive P2 (eq 
(5)) from P?,; and Py (eq (3) and (4)), although the 
individu: al b ge tone a P, and Py, do refer to 
processes in channel I and in channel IT that are 
comeinely vendieta. 


3. Rigorous Derivation 


In determining the dead-time loss and accidental- 
coincidence corrections that must be applied in the 
method of 42-crystal-scintillation counting, it 1s 
necessary to apply a rigorous analysis of the events 
occurring in the two multiplie r-phototube channels 
that lead to the coincidence analyzer. This Is neces- 
sitated by the fact that a number of the events 
occurring in one channel are linked to those occurring 
in the other, and it may be therefore no longer valid 
to apply the kind of considerations that would be 
possible if all the events were completely randomly 
distributed. 

Let us assume that NV, counts per second, the 
so-called true counts, arise in the composite crystal 
scintillator due to the radioactive source and _ the 
background. Then, of these counts, a number Ny 
per second will be recorded by the coincidence 
analyzer. Further, let Ni and NZ be the number « 
noise counts per second passing along channel I and 
channel II, respectively, and_ let N.. be total 
number of coincidences recorded per second | Vv the 
coincidence analyzer. 
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also necessary to consider, in each 
be conveniently designated as 
which we will assume to be 
per second. These erise in 
the following manner. A true count, if it can pass 
through both channels unimpeded, arrives at the 
coincidence analyzer as a ‘“‘perfeet’’ coincidence with 
no delay whatsoever and, as an “N,,”’ event, will be 
recorded. If, however, the pulse due to a true 
count happe ns to coincide with the gate of a preced- 
ing noise pulse in one channel, the true count. in 
that channel will be lost and the pulse in the other 
channel due to the true count will proceed to the 
coincidence analyzer linked to the noise pulse that 
has “consumed” its partner. A divorced true count 
is thus always linked temporally to a noise count 
in the other channel, and such divoreed true counts 
are therefore no longer randomly distributed with 
respect to the other counts. 

The diegrams in figure 1 will help to illustrate 
the relationships of the various kinds of counts one 
to another. In these diagrams the pulses of noise 
counts are shown above the time axis of each channel, 
whereas the pulses due to true counts are shown 
below the axes, although in reality they are @X- 
perimentally indistinguishable. Figure (a) illus- 
trates the circumstances in which a divorced true 
count (V; per second) can arise in chennel I, figure 
l(b) the circumstances giving rise to a divoreed 
true count (Nf per second) in chennel IT, and figure 
l(c) those in which the true counts ean be lost in both 
chennels I end Il. In all eases the length of the 
pulse is chosen to signify the time 7, in which the 
electronic gates in either channels I or II will be 
closed to the transmission of further pulses. 

What then are the numbers of divorced true counts, 
Ni and NY per in chennels IT and II, 
respectively? 

These correspond to the numbers of overlaps be- 
tween true-count pulses and preceding noise-count 
pulses in the other channel, for times greater than 7¢, 
the resolving time of the coincidence analyzer. 


It is, however, 
channel, what may 
“divorced true counts” 
respectively Ni and NY 


se ‘ond, 


$0400) 


ne 


If a pair of true counts is impressed on the cir- 
cuit within a time interval less than 7, after a single 
noise pulse or a coincident pair of noise pulses, the 
coincidence analyzer will accept these as a true co- 
incidence. The effective noise-pulse dead time for 
the cre nae of divorced counts in the other circuit 
is thus ry, 

In shaninel II, the noise counts Nf per second 
give rise, therefore, to an effective total dead time 
of Ni(7g—rT) per second for the creation of divorced 
true counts (V{ per second) in channel I. During 
this time N, true counts will be submitted to chan- 
nels I and II. Therefore, if there were no inter- 
action between channels I and IT, 


divoreed true counts in channel | 


(rg—Te) NAN; per second. (16) 
Similarly, 
divoreed true counts in channel II 
ayaa a 
(te—Te) NaN; per second: (17) 


These expressions as they stand give, however, 
somewhat higher values than N{ and N{, because 
they also include cases of overlap, as in figure 1(c), 
in which both true-count pulses are rejected by 
two preceding noise counts. 

The numbers of overlaps between a pair of true- 
counts and a preceding noise count in channels [ 
and Il are, respectively, ro )NAN, and (tg—T¢) 
Ni; per second. 


The time in which a second noise count in channel 


(T, 


Il (fig. 2) must occur to eliminate N/’ varies therefore 
from 0 (fig. 2(a)) to (7,—7.) (fig. 2(b)), giving an 
average time per overlap of 3(tg—T-e). 


Thus, considering figure 2, the number of such 
divorced true counts passing along channel II is 
(te—Te) NaN: per second (eq 17), and the average 
dead time for each of these events for the acceptance 
. preceding noise count in channel IT which will 
pulse is 3(t,—7,-). This corre- 
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FIGURE 2 Overlap helweenl 


sponds for such events to a dead time per second of | becomes available for acceptance of 


B(te—Te)?NiN,. With Ni noise counts per second 
in channel I], the number of triple events causing the 
loss of both of a pair of true counts in channel I and 
channel I] a(te—Te) NiNiN; per second. It 
is seen that the same result will be obtained by 
proceeding via ‘“N,”’ events in channel I (eq 16) and 
the Ni counts in channel I. 

We thus derive the following corrected values for 
the divorced true counts in channels I and IT: 


Is 


Ni=(te— te) NENG {1 —43(tg—7-) Nz} per second, 
(1S) 
and 
Ni=(ty— te) NN: {1 —4(7g—7-) Nf} per second. 
(19) 


It is now possible to proceed to the computation 
of the true counts lost as a result of the events that 
are summarized in figure 1. In the first place, every 
“Ni” and “NY” event corresponds to a lost true- 
count pair. In addition, true counts can be lost in 
the dead times initiated in each channel by the 


divorced true counts, Nj and NY per second. The 
: . . CATS SS cc ATM ss 

dead times involved in N,” or “N, events vary 

from 7, to r,—7T,, as illustrated for the “N?”’ events 


in figure 3. 

The dead time is only that which swueceeds an 
‘“N.”’ or ee pulse. In no case does the preceding 
noise count contribute to the dead time, for, from the 
very nature of the divorced true counts, an earlier 
true-count pair following a noise count merely creates 
an earlier divorced true count. 

It is a curious fact that, as shown in figure 3(¢), an 


intermediate time interval, varying from 0 to ze, 
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ie-countl p tine 














and p eceding noise-count pu 


true counts 
immediately following the noise-count dead time in 
channel I. This is effeeted when the channel-I pulse 
of the true-count pair registers as a coincidence with 
the “nN” pulse at the coincidence analyzer. The 
fraction of the ““NU’ or “NY” events in which this 
window becomes available is 7,/7,, and its width is on 
the average equal to $7... Thus, the true counts that 
can be accepted in this window for “Nj and “NY?” 
events are, respectively, A 7? reN.N, and 472 tN N, 
per second. 

It therefore follows that the total true counts lost 
in channels I and I] due to the dead times succeeding 


° 
. 


oN and “NY” events are respectively NUN,(7,.- 
, ad , , 
tr? T,) and N.N.(t.; tr? Tg) per second. 
Pa , a - 
In addition to losses of true counts N,, N,, 


NUN. (rg—422/7,), and NYN,(7,—3 true counts 
can also be lost in the triple events illustrated by 
figure I(c These number 4(r,—7¢)?NANIZN;, per 
second. These triple events in which a pair of true 
counts is completely eliminated by fo noise counts, 
one in each channel, are illustrated in greater detail 
in figure 4. 

As can be seen from inspection of figure 4, the 
dead time due to the second noise count varies from 
r, (fig. 4(a)) to the average of 7, and 7,, or 3(tg+Te 
(fig. 4(b This is equivalent to a total average 
dead time of 3(37,+7,). The true counts lost in 
triple events such as that of figure 1(c) are therefore 
g(Te—Te)” (S7_- te) NNN? per second. Therefore, 
the total true counts lost per second as a result of 


3Tc/Ts), 


Im the tual experiment the resolving time f the amplifier was eater than 
Tr so that coincidences such as are indicated in figure 3(« ao not occur Figure 
3(c) is vacid only for very sharp pulses whose duration is much less than + 
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losses in both channels I and IL due to overlapping | > Nu 4d Nt —3 )N") 
true-count and noise-count pulses are | aay, Ste Te ali — 3\ Te Te)2Vn 
rNu N.+N, N.( 14 2 ) (MAND) +N, (] 4( tT, OND} {IEMs 1 ‘ )} 
Ta Tg 
1 2 "2 «p j ye ” r2 vey vr , r ‘ 
s\T¢ Te) \OT, Te \ nN tte T 1 ( T. Te) (OTe t.)N,N,N:- (21) 
© SNucx (N’ +N’ \ 1 Te N,. is the number of true-count coincidences per 
eal = 4 4 aVr | avr Tg 2 | . . . 
Ve | second that are registered by the coincidence analyzer 
| and therefore pass through channels I and IL un- 
+$(7,—Te)? (3724 r.)NINZN?2. (20) | impeded. Each of these counts imposes a gate 
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dead time of 7,, and gives rise therefore to V,.7,.N; 
losses per second. Therefore, the total of all true 
counts lost from any cause whatsoever is equal 
a ta: 
Nie ite T S WN a. 
ry ’ . 
The true counts per second, N,, are therefore 
given by 


i 
, Nee 
N, a 
\ Nu (22) 
l cel gn nF 
N, 


Ni. being the number of true-count pairs that pass 
through channels I and Il unimpeded. \, is, how- 
ever, the fotal number of coincidence counts re- 
corded per second. Thus the difference between N,, 
and N, will be essentially due to accidental noise 
coincidences that are of the order 27.0Vi—N? 
(N{—N/,) per second. Thus, 
Nie~N.—217-(N,—N,)(Ni—N,). (: 


n 


~ 


The reason for subtracting NY from NY and N/ 
from Nf in the aecidental noise-coincidence term 
is that any noise count associated with an “.NV{" er 
“NT” event can never contribute to en aecidentel 
coincidence. This xecidental coincidence correction 
in eq (23) is, however, en overcorrection in that true- 
count pulses can arrive during the time the accidental 
noise coincidences are bet Ing accepte “dl by the coinei- 
dence analyzer. The time during which such true 
counts can be aecepted as coincidences is seen by 
inspection of figure 5 to vary from 0 to 2r,, giving 
an average of r,. 

The number of coincidences between accidental 
noise coincidences and true-count pulses is thus 
equal to 2r7.(Ni, Ni)(Nz N,)7N, per second. — In 


other words, the accidental noise-coincidence term 


must be multiplied by a factor equal to (1—7.N,) to 
allow for true counts coinciding with accidental noise 
comcidences. Thus, as a closer approximation, 


Nie~Ne—2r(Nan—Ni)(Na—Ni)(1—reN,). (24 


There is, however, still one other possible source of 
coincidence counts that cannot be overlooked. If we 
consider figure 1(a) in detail it will be seen, as shown 
in figure 6, that a window also exists for a very small 
number of coincidences of “N/"’ and “NY events 
With noise counts, just as a window was previously 
noted to exist (fig. 3(¢)) for the acceptance of true 
counts. 

The numbers of “Ni" and “NN?” events that 
een be available for such coincidences are equal 
respectively to [r.(rg—re)INi and [re (te— red |NG 
for 27.<7, or Ny and N{ for 27, The average 
time that the window is open for the acceptance of 
a noise count in coincidence is }7,., provided that 
2re<ity, OY 3(3r—e—T7,) if 27.>7,. The numbers of 


g- 


T he 


/ 


such coincidences for “Ni” and a? events, 
a ging. ere therefore }$[r2 (7, JINN and 
+ [72 Ty IN: NN, n per second as, in re ne ral, T o< Tyg. 
As these coinek te ‘neces can also be in coincidence with 
true-count pulses, they too represent an overcorrec- 
tion and must be multiplied by the same correction 
factor (1 ToN;) used to correct the accidental 
noise-comeidence counts of eq (24 
We can therefore now write exactly that 


Nie=Ne— <{ 2t(Na—Ni)(Na—N)) 


| T. ee —— - 
tu (NL.NAt+N,.N,) ] r. (25 
6 te Ti 
Here again, however, because the resolving time of the amplifier w reater 
than re, these coincidences could not, in practice, oecur 
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FIGURE 6 “Window” 


The true counts per second, N,, are now given by 


N, ENG - 


t 


; Nu 
where .V,, is given by eq (25) and ~ bv eq (21). 


N, 

The expression for N, can therefore be written 
down exactly for those who desire to experience such 
an exercise or whose experimental conditions demand 
it. In the interest of economy of space, however, 
let us at this point make an approximation that is 
justified by the physical quantities involved under 
the experimental conditions [1], namely, to neglect 
all terms involving higher powers than 7Z (7,.~10 
usec) and +r, (about 4 in eq (21). Then, by 
substitution from eq (21) into eq (22), 


usec ) 


Nt 
1—N et Te—Tte IN N, Te—Te)N, +N ites} 


26 


N,$t] 


It is seen from eq (18 and (19) that substitution 
for Vj and N7 in eq (25) for N,. will introduce terms 
involving +r,.(7,—T,.) and r.)*?, which can be 
neglected by comparison with 27,.N,N%. Equation 
(26) can therefore be further simplified to give 


TelTe 


N.—2reN,N 
[IN 


Ws l T (N, 2r.N N (T¢ (Tr T. iN N 


T N, s 


Te)} 
As a further simplification, because 
NitNimN(D + ND) —2(N.—27eN Nn), 
if all products of 7's contaiming T. are neglected, 
N.—2reNaN 
+N 


"1—7,N.—(t.—7. {NO 11)—2N.}+72N,N2 


(28) 





TIME—> 











for coincidences between noise counts and divorced true counts, 


In order to determine Ny, and Nf a blank experi- 
ment can be carried out with the crystal removed and 
with an opaque obstruction between the two multi- 
plier photo tubes, to prevent phosphorescence in one 
photocathode causing a simultaneous signal in both. 
Then .\, is zero and the noise counts will be registered 
directly by the sealers in each of channels I and II. 
By inserting the crystal and carrying out experi- 
ments with and without a source, the values for N,, 
for both the source plus background and_ back- 
ground, respectively, can then be determined in 
terms of the number of coincidences and the known 
quantities r, and 7, and the noise counts. 

It is interesting to note from an inspection of eq 
(18) and (19) that if r,—7,, the numbers of divorced 
counts N{ and NY become zero. For such a condi- 
tion, S.N,, (from eq (21)) also becomes zero and eq 
(22) gives without any approximation 


N.—2teNaNa(l —teN2) 


N, ; aE NT =—~e 29) 
1— Nore+27r-NaNa(1 —7reN;) 


In considering the case for 7z=7.¢, it is also necessary 
to use the form of eq (25), for Ni., where 27, >7¢, 
namely, 


Ni=Ne— {2r, CN, 


= —— 
Tare” tz (NINat 


Na) (Na— Ni) 


NIN,) (1—r.M.) (30) 


in which, however, the second term in the first 
bracket becomes zero when rz=7, by virtue of both 
Ni and Nf becoming zero. However, this is never 
experimentally feasible, for (a) if 7, is inereased to 
equal 7,, the accidental rate will be too high, and (b) 
if r, is decreased to equal 7,, the non-randomness 
tends to disappear and the rigorous derivation be- 
comes unnecessary. 
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Limits of error are indicated on only two points, 


4. Comparison of Rigorous and Approxi- 
mate Derivations 


It is of interest to compare experimentally the 
validity of the corrections required by the approxi- 
mate derivation, eq (12), and the rigorous deriva 
tion, eq (28). The arrangement shown schemati- 
cally in figure 7 was designed to test directly both 
derivations over the range of phototube noise en- 
countered in the experiments. 

A sandwich source of TP, covered by a hemi- 
spherical aluminum reflector, wes counted by photo- 
tube S. These counts, and the counts from 
phototube S, were fed into channels I and IL simu/- 
taneously, thus providing a source of “true” pulses. 
Separate phototubes, N,; and N»., were then used as 
independent sources of random noise, and pulses 
from N, and Ne were fed at the same time into 
channels I and II, respectively. The counting rate 
from S was maintained constant. 

The voltages across the dynodes of N,; and N» were 
varied stepwise so as to increase the random noise 
from 0 to 11,000 counts per second, covering the 
most extreme noise measurements encountered in 
the actual experiments. The noise pulses from N, 
and the ‘true’? pulses from S passed through channel 
I, and similarly the noise pulses from N, and the same 


HOUSE 


“true” pulses from S passed through channel IL. 
N(D, NCI), and N,, were then analogous to the 
values obtained in the actual experiment, except in 
this case V, was known exactly (by counting S alone) 
and Ny and Nf were known exactly (by counting N, 
and N, separately). In figure 8 are shown the values 
of N, (approximate) and NV, (rigorous) obtained by 
means of eq (12) and (28), respectively, as functions of 
N(1)+N(11), the sum of the single-channel counting 
rates. As can be seen, both methods of correction 
are applicable over the range of phototube noise en- 
countered in the experiment. 

The foregoing derivations will be applicable to 
any system where coincidence between two counters 
is used to select particular events from a_high- 
background environment. 


We gratefully acknowledge the assistance of Carter 
C. Smith in carrying out the experimental comparison 
of the rigorous and approximate derivations. 
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Entropy Changes in Rarefaction Waves’ 
Robert F. Dressler 


Frictional flow of a polytropic gas is investigated for centered rarefaction waves. En- 
ergy balance is maintained by assuming that mechanical energy loss due to frictional force 
reappears as heat. Entropy behavior at the wave front is discussed. The nonhomogeneous 
linear system of partial differential equations with variable coefficients for first-order fric- 
tional effects is derived. Utilizing the geometrical similarity of the Mach lines, the fune- 
tional form of these quantities is ascertained, which permits explicit solutions of the boundary 
value problem. The expansion procedure is singular at the wave front, and results are not 
applicable there. First-order effects are expressed as polynomials plus terms singular at 
the wave front. Results are compared with expressions obtained when heat generation 
due to frictional force is neglected. 


1. Introduction 


Equations for one-dimensional, unsteady flow of a compressible gas permit explicit solu- 
tion for centered rarefaction waves, such as occur in shock tubes or pneumatic control circuits. 
For some applications, however, it becomes important to include also the complications re- 
sulting from frictional dissipation. This matter is discussed here in two ways, first with the 
simplifying assumption that the frictional effects consist merely in a retarding force, and 
second, by the more complete consideration that includes the resulting heat generation and 
accompanying change in entropy. This means that the mechanical energy destroyed by the 
retardation force reappears as heat energy, resulting in additional expansion. The simpler 
flow description without heat generation will be referred to as model A, and the second, more 
accurate, one as model B. The equations for steady frictional flow in model B have been 
studied and yield energy theorems and an extended form of Bernoulli’s theorem [1],? but knowl- 
edge of the unsteady flows appears to be restricted mainly to a few isolated numerical calcu- 
lations. 

A centered rarefaction wave is perhaps the most basic pattern in unsteady flow; the present 
paper derives solutions for first-order dissipative effects in this unsteady wave using the full 
system of equations of model B, and compares results with corresponding ones for model A 
(which were previously published by the author [2] as a supplementary result in a study of 
shock tubes with varying cross section). 

A recent paper by Ludford and Martin [3] has discussed anisentropic effects in centered _ 
simple waves. They consider flows where the specific entropy varies from particle to particle, 
but remains constant for a given particle. In frictional flow, however, the specific entropy 
for a given particle must vary with time, and a centered rarefaction wave cannot remain @ 
simple wave. 

Consider a polytropie gas at rest for time ?<0 in a duct of uniform cross section, filling 
the duet to the left of a diaphragm where 7<0, and with a vacuum initially at 7>0. Heat 
exchange between the gas and the duct, heat conduction and radiation, and viscous effects are 
disregarded. When the diaphragm is removed at f=0, a centered rarefaction wave begins. 
Dimensional quantities are denoted by bars in order to define the unbarred dimensionless 
quantities to be used subsequently. Before flow, the gas state is described by sound speed 
@, mass density po, and temperature 7p; after flow we have velocity 7Z, local sound speed e¢, 
mass density p, and pressure Pp. 

For model A, assuming no internal heat generation and no entropy changes, the equation 
of state can be written in forms 


p=A%", S=yA*p'=7RT, (1) 


his research was supported by the Office of Naval Research, USN 
Figures in brackets indicate the literature references at the end of this paper 
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where y is the adiabatic exponent, R the gas constant, 7 the temperature, and A* is a known 
constant (consistent with the assumption here of no entropy change). In this simpler model, 
the friction produces merely a retarding force per unit mass of magnitude 7 where X is the 
known friction factor. There is some question concerning the constancy of this coefficient for 
frictional flows that are highly unsteady; this question is discussed by Schultz-Grunow [4], 
Jenny [5], and others. Here it is considered to be a constant. 

Introducing dimensionless (unbarred) variables defined by 2=7/D), t=(@,/D)t, u=W%, 
C=T/E, p=p/ po. X=DX where D is some length associated with duct cross section, the momen- 
tum and continuity equations for model A are 


~ 


u,+uu,+—— ce 
TI 


for the two unknowns wu (z,f) and ¢ (r,t); and the energy equation is neglected. 

In deriving equations for model B one must take account of the mechanical energy lost 
through the action of the retarding force \wju. The action of a body force cannot itself di- 
rectly change entropy; however, if a term is included in the energy equation for an amount 
of generated heat equal to the lost mechanical energy, this heat will change the specific en- 
tropy. In this way, the energy equation becomes expressible directly in terms of the friction 
coefficient. The equation of state for model B now includes the entropy dependence, 


p= A(s)p", A(8) = (¥—1)e?~ 0", e=yRT(p,3) (: 


i 


~ 


where ¢, is specific heat at constant volume, § is the specific entropy (entropy per unit mass), 
and 8, is a constant with a value depending unon the entropy level adopted. Letting dimen- 
sionless specific entropy s be defined by s=(7)/@2)%, the momentum, continuity, and energy 
equations for model B become, respectively, 





9 ) 
u,+ uu,+— ce, hu u +c?*s, 
Y 
2 2 u 
cu,+ uc.+ ( AY } 
Y l 7 l ( 
ul 
x WS X 
a 
4 
for the three unknowns u(r,f), e(7,f), and s(t). We set s=s=0 for the gas when initially at 
rest, and the equality 
\* 1(0) \4 le ° (5) 


determines the constant 8». 
2. Characteristic Equations 


The Mach Jines for mode! B are defined by the system of characteristic equations equivalent 
to the totally hyperbolic svstem (4). These consist of three families of real curves with direc- 
tions given by 





dy 5 
u—cC 
( di ( 
ly 
C9: a u-+c > (6) 
dr 
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which are respectively the backward sound path, the forward sound path, and the particle 





trajectory (see fig. 1). The corresponding differential relations along each family are 
2 a 
te d( a ;¢ )te ds hulu ( 1+ fy—1] - ) dt 
) 
(2 d( sia” - | c )—e ds AU U ( ] ly 1] ~) at > (7) 
( ds=h - dt 
Cc a 


These equations are used in [6] as the basis for a discussion of a numerical computation on a 


rarefaction wave in a shock tube. 


3. Entropy Behavior at the Wave Front 


From the above equations and the equation of state (3), the behavior of s at the forward 
wave front where p=0 can be inferred, consistent with the assumed model for the flow. For 
notational simplicity and to permit numerical results, we now specialize the gas to be air with 
y=7/5; however, the same arguments will apply and analogous solutions can be derived for 
any other admissible value of y. 

If % is chosen as indicated and the relation y=1+ 7? ,, is used, we obtain the dimensionless 


ie .v=i 


equation of state in the form c?=¢«7°~'*p’~', which for air is 
gag" . (8) 


On the wave front curve where p=0, this permits two possibilities: either (a), s remains finite 
and hence e=0 there, or (b), s becomes infinite and ¢ is not known in advance. — If s were finite 
and continuous and ¢ vanished there, then by the third equation in (7), the rate of change of s 
along the wave front would be everywhere infinite as «#0 there. This follows from the fact 
that the wave-front locus is also a particle trajectory and therefore a &3 Mach line. Therefore 
< must actually be infinite on the wave-front curve, and ¢ is not determined there by this argu- 
ment. Because s is entropy per mass, this does not necessarily imply that entropy per volume 
becomes infinite at the wave front for this model. 


4. First-Order Effects 
By setting up asymptotic expansions of form 
r(r,t,rA)~r%(r,t)+ VA(7,OA+ . . . (model A) 


(8) 
r(a,t.A)~r(r,t)+ V%(r,trA+ ... (model B) 
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a perturbation scheme will be defined to give information for small values of A. The quantities 
v(z,t,A) represent each unknown w,c,s, and their various first derivatives. The unperturbed 
centered rarefaction wave is the known simple-wave solution 


u®=5(1+2/t) 
c’=34(5—2/t) (10) 
g’==§ 


which when A=0 satisfies system (2) for model A and system (4) for model B. The first-order 
terms for model A are 1, C*, (and S*=0) and have been published in [2]. Now we proceed 
toJsolve the more complicated problem for (®, C®, and S® and to compare these with the 
previous results. 

The expansion s~0+S®(a2,t)A+ . . . shows at once that the perturbation procedure for 
model B cannot be expected to retain validity in the neighborhood of the forward wave front, 
since s(z,t,A)= © for any A>0 on the wave front whereas s(7,t,0)=0 there. The results 
therefore will not apply in the neighborhood of the forward wave front where a type of boundary 
layer effect occurs. A separate analysis would be required in that region, possibly similar to 
the application of the Pohlhausen method as used by Whitham [7] for a water wave problem. 
There would be little justification for it in this problem, however, because the model assumed 
for frictional dissipation would in any case become inaccurate in this region of extremely low 
density. (A singular situation does not occur in the perturbation for model A. In that case 
there is no edge layer effect, and results hold throughout the full wave; but on the other hand, 
the basic model itself is everywhere less accurate than model B.) 

When expansions of type (9) are inserted into system (4), after using (10), the first-order 
terms satisfy the system: 


* 
(5t—z)* 30(¢+ 7) ,- : 30(5¢—2z) / 30 25(t-+ 2) 
S8— U7®—36/* cPr—— ['® Ce4— 0 
r? : f Mies f ce" © 2 
(5¢—z),- Q2A(t+.27) 5 W75(t+.s) 
(B+ CB+30CR——_ ( 84+— (®— 0) > (11) 
t : f a t 617(5t—.z) 
d( t+ 7) SB4 g.gp—12ol! t r) 0. 
f t(o5t—z)- a 





The solution of these equations that is sought must satisfy the conditions /®’=C®=S®=0 on 
the straight line t=—z, (t>0). This nonhomogeneous linear system has the same character- 
istic curves as the solution uw’, ce’ plus its particle trajectories. By studying these curves one 
can infer an easy procedure for obtaining the desired functions. These curves in the flow 


wedge are 


: rz=mt, —loimss5 > 





where each equation defines a one-parameter family of curves over the range of a parameter as 
indicated. The ¢$ and ¢f curves emanate from the back wave line s tat a point (—a,a) 
for any a>0 (see fig. 2). Writing ¢) as (r/a)=m/(t/a), one sees that the coordinates of an 


intersection point of a fixed ¢) curve with any 


Therefore as the parameter @ is varied, the quantities t/a and s/a at such intersection points 


curve must be proportional to the parameter a. 


remain constant as we progress along a fixed ¢f ray. Likewise an analogous situation will 


268 
































FIGuRE 2. 








apply for intersections of ¢ curves with a fixed ¢? ray. This geometric similarity of the char- 
acteristics can be employed to find the functional form of the desired solutions (7%, C®, S®. 
The characteristic relations associated with (12) are 


5t—z 25(2-+-ty 25(2r-+-t)* 5 U®—c® a 
re ({ B 5 B ) N\B — === 
a: d ' + Gt )d ( 18¢2(24—5t) 360? )at Ss 4 dt 
5t—2zr 25(a+t)* 25(a+ t)? 
oO. ET B45078)—(” 1 f (13 
- * ok. ( Gt )e ( i8fr—t) 3682 )d F 
25(t+ 27) 
Yo: dS® : “ v) at. 
sy 6? (5t—az)- . 





The last equation does not contain unknowns on the right side and could therefore be integrated 
along any ¢3 curve. If this were done, using the third equation in (12) and the geometric 
similarity, it is easily seen that it would lead to a relation of the form S®/a=F;(t/a) when 
S®—0 at the origin. Because t/a is constant along a ¢ ray, S® is proportional to a, and hence 
to ¢ there; the functional form must then be S®(z, t)=L(m)t. Using this result for dS® in the 
second equation of (13), by analogous argument applied to the ¢ curves, it follows that (78+ 
5C®—F,(m)t. Finally, we see that if we integrate the first equation of (13) along any ¢ ray 
(where m=constant), we would have ((/®—5C®) proportional to ¢ if the quantity 5(U/®—C®) /3t 
were a constant along the ray. Then it would follow that (’® and C® would separately be pro-, 
portional to f also. But in that case, 5(17®—C*)/3t would actually be constant as desired. 
There these considerations indicate the existence of a solution to (13) and (11), vanishing at the 
origin as required, with the functional form f(m)t for each unknown Ll’, C®, and S®. This 
information now permits immediate solution of the original system (11), as follows: 

First, system (11) is transformed through use of the new independent variables m and ¢, 
putting / (2, t)=U7*(m, t), and likewise for the other unknowns. Then setting l’*=A(m)t, 
OC*—k(m)t, and S*=L(m)t leads immediately to a nonhomogeneous system of three ordinary 
differential equations for h, k, and LZ. Although this system contains variable coefficients, it 
can nevertheless be solved exactly. The particular solution satisfying the conditions h(—1)= 
k(—1)=L(—1)=0 is then finally obtained to give the first-order effects for model B. 


U® 2517 35 ; 180 7 

- (5—m)y—7(5—m)*+ 54(5— m)—172+- 
a T6—m) 
C® 5T(5—m) , 35.. — : - , 1868 90 71i+m)y 

— he nn (5—my—15(5— m)*+108(5— m)+(1+my—— — — > (14) 

t 48] 72,576 »4 i (5—m) (5—m) 
S® (5—m)° l 108 27 

2125 ——(5—m)+3-—=: - - | 
kan Ree 5 a m)* (6 -m)° . 
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For comparison of these results with the corresponding expressions for the simpler flow 


model, we list below the expressions obtained in [2] for model A: 


(vA 2512 5 
(5 m) AS m)-+36(5 m) S4 

t 144.9 

c= 5 . 2 7 as ’ “ os 
; 144 g m) 10(5— m+ 72(5— m)— 168+ 3(m+-1) P (15) 
S4 

0. 
f A 





The expressions in (14) have singularities at m—=5, at the forward wave froat for the uaperturbed 
flow. This is to be expected in view of the prey ious remarks conceraing the entropy behavior 
at the forward wave front that makes the perturbation procedure singular. On the other 
hand, the expansion for model A does not possess this singular property, and the expressions 
(15) do not become infinite at m=5. A comparison of (14) with (15) must be restricted to a 
region that excludes the forward wave-front zone. Figure 3 gives the numerical evaluation of 
the quantity S®/t. When drawn on a semilogarithmic plot, this curve is almost linear over the 
middle range 1<m<4, and for convenience can be accurately approximated there by the 
expression .63 exp (1.95 m). Figure 4 compares the behavior of ('®/t with (4/f. In the first 
half of the flow zone, the two models predict almost identical results, but further on the velocity 
becomes higher for model B. In equations (4), as s, >0 for this case, the effect of the entropy 
term in the momentum equation is to diminish the tota! resistive force, and this by itself would 
increase the velocity. In the continuity equation, however, the extra positive term on the 
right side would have the effect by itself of decreasing the velocity in the back portion of the 
flow zone, because this term is analogous to the term that would be present for a converging 
duct in quasi-one-dimensional flow, as derived in [2]. Evidently the fact that the velocity 
term for model B is slightly below the model A curve in the beginning of the flow zone, is 
due to the extra term in the continuity equation dominating over the entropy term in the 
momentum equation until s, grows sufficiently large to raise the velocity. 
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The first-order corrections for the local sound speed ¢ are shown in figure 5 for the two 
models. In this case there is no resemblance in the behavior of the two curves, as the sound 
speed, Which is proportional to the square root of the temperature, is more directly sensitive 
than the velocity to the internal heat generation assumed in model B. 

Finally we use our results for sound speed and entropy in the two models to calculate the 


corrections for the dimensionless density p= p/p). The dimensionless equation of state of air 
for model A corresponding to (1) is p=c®. After inserting expansions here, the coefficient of 


in the expansion of type (9) for p is then equal to 5(e°)'C*. This quantity, after division by ¢, 
is shown in figure 6. For model B, when corresponding expansions are inserted in (8), the 
coefficient of X in the p expansion is equal to 5(e°)* [(C®—(7e°S®/25)]. This quantity, after 
division by ¢, is likewise shown in figure 6, and one observes that the widely different results 
of the two models for sound speed and entropy nevertheless combine to produce density cor- 
reciions that are in close agreement. The pressure and temperature behavior for the wave can 


be obtained directly from (3) using the results already presented. 
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Statistical Investigation of the Fatigue Life of 
Deep-Groove Ball Bearings 
J. Lieblein and M. Zelen 


Fatigue is an important factor in determining the service life of ball bearings. Bearing 
manufacturers are therefore constantly engaged in fatigue-testing operations in order to 
obtain information relating fatigue life to load and other factors. Several of the larger manu- 
facturers have recently pooled their test data in a cooperative effort to set up uniform and 
standardized ball-bearing application formulas, which would benefit the many users of anti- 
friction bearings. These data were compiled by the American Standards Association, which 
subsequently requested that the National Bureau of Standards perform the necessary 
analyses. This paper summarizes the principal results of the analyses undertaken by the 
Bureau, and describes the statistical procedures used in the investigation. 


1. Introduction 


1.1. Statement of Problem 


The experience of ball-bearing manufacturers over many vears has led to the acceptance 


of an equation of the form [15, p. 15, eq (53)] ! 
L=(C/P)?, (1) 


relating fatigue life Z to load ? when other factors are kept constant. In the above equation, 
(‘is termed the “basie (dynamic) capacity,” and is defined [15, p. 48] as the constant bearing 
load (in pounds) that 90 percent of a group of similar bearings can endure for one million 
revolutions under the given running conditions. 

The quantity Cin eq (1) depends upon the characteristics of the bearing type, as indicated 
in [15, p. 32, eq (120)]. When the expression cited is substituted in eq (1), the fatigue-life 
formula for ball bearings takes the form 


I fZ2°%DdD2(4 COS a )%3 a 9) 
y P (2 


The svmbols are defined as follows: 

Z—number of balls. 

1),=ball diameter in inches. 

/=number of rows. 

a=contact angle. 

P bearing load in pounds. 

L—number of million revolutions that a specified percentage of bearings will fail to sur- 
vive on account of fatigue causes. If the percentage is 10, then L=Lyp, and is 
termed the rating life; if the percentage is 50, then L=L,, the median life. 

p, @), @2, Ay, fe ave taken as unknown parameters whose values have to be estimated from 

given data. 
Since 7=1 and a=0° for deep-groove ball bearings, with which this paper is exclusively 
concerned, the life equation that will henceforth be considered takes the form 


L E => | (2a) 


Figures in brackets indicate the literature references at the end of this paper 
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The main goal of this investigation was to determine the “best values’ of the unknown 
parameters in the life equation from the experimental data. One of the major problems was 
to determine the value of the exponent p, as there was disagreement within the ball-bearing 
industry whether an appropriate value for p was 3, 4, or some other value. 


1.2. Description of Data 


The data available for analysis consisted of sets of records summarizing endurance tests 
for deep-groove ball bearings. These tests were carried out over a period of vears by four 
major ball-bearing companies. In the interest of trade anonymity, these companies will 
henceforth be designated by A, B, C, and D. Each endurance test consisted of a number of 
bearings of the same type (the number varying from test to test), which were tested simul- 
taneously under the same load and running conditions. Table 1 summarizes the number of 
test groups of data for each company. ‘The data from company B were sufficiently extensive 
to permit a further breakdown into three bearing types, here denoted by B-1, B-2, and B-3. 


TABLE | Summary of ball-bearing data 
Number of test Total number of 
Companys groups bearings in test 
groups 
50 1, 259 
B 148 3, 289 
Type B-1 37 
Type B-2 Of 
Type B-3 17 
C 12 201 
DD ; 109 
Total (all companies 213 1 O48 


The worksheets, summarizing the tests, recorded the number of millions of revolutions 
reached by each bearing in the test group before fatigue failure. Information was also given 
for those tests terminated before all bearings in the test group failed. In addition to the test 
results, the worksheets included information on the characteristics of the bearing type (e. g., 
values for Z, D),, 7, a) and load 7, as well as other items of descriptive and identifving informa- 
tion. A specimen worksheet is reproduced in appendix A. 

All necessary quantities for evaluating the unknown parameters in the life equation (2) 
were given directly on the worksheets except the fatigue life 1.2) This quantity can be estimated 
from the observed fatigue lives of individual bearings within a test group. As already noted, 
two concepts of fatigue life are used for L, namely, the rating life Ly, and the median life Leo. 
Separate analyses have been carried out with regard to each throughout. 

Appendix A summarizes the data taken from the original worksheets that were used in 
the statistical analysis. Also given are the computed values for Ly, Ls, and the “Weibull 
slope” e (which relates to the dispersion of fatigue lives). The methods for obtaining these 


quantities from the bearing data are given in detail in appendix B. 


1.3. Assumptions for the Statistical Analyses 


All conclusions reached in this report, and all statistical analyses employed, are based upon 
the following principal assumptions: 

(a) The life formula (2) is the proper functional form for describing fatigue life in ball 
bearings. 

(b) Differences in the measured life of bearings classed as identical, tested at the same load, 
reflect only the inherent variability of fatigue life, and are free from systematic errors that may 
arise from different test conditions, materials, manufacturing methods, ete. 


Certein estimates for / ind J had been entered on the worksheets for many of the test However, these were not regarded 
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(c) All the bearings in a test group can be regarded as a random sample from a homogeneous 
population of ball bearings. 

(d) The probability distribution of the number of revolutions to fatigue failure is of the 
same form for each test group, although its parameters may differ from group to group. 

(e) This fatigue-life distribution is of the type known as the “Weibull distribution.” 

The purely statistical assumptions, (c) to (e), served as the basis for the determination 
of Lio, Ls, and e for each test group. Assumption (e), however, is not involved in the methods 
used to evaluate the parameters in the life formula (2) from given values of Lyo or Ls. A 
different assumed form fer the distribution of fatigue life might give somewhat different values 
for Ly) and Lj», but the same methods could then be used to evaluate the unknown parameters 
in the life formula (2). 

Other assumptions of a more technical nature were necessary in the course of the analyses. 
These are discussed in appendixes B and C. 

As in all cases where inferences are made from given data, the conclusions reached here 
pertain only to the population from which the given data can be regarded as constituting a 
random sample. 


2. Outline of Statistical Analyses 


The statistical analvses were divided into two phases. The first phase considered the 
problem of finding estimates of Ly, Leo, and the Weibull slope e from the given test data; the 
second phase used these estimates of Lj). and L;, to evaluate the unknown values of the param- 
eters in the life formula. 

2.1. Estimation of Lj») and Ly» 


The quantity 1 depends upon the existence of an underlying probability distribution of 
bearing lives. Selection of a distribution or population is equivalent to specifying the proba- 
bility that a bearing selected at random from such a population will survive any given number 
of revolutions, L, or, conversely, that if ¢ is a specified probability, then L is the life period 


- 


that will be survived with this probability, e. g., 


.90 for L=TL yp. 
Probability {life>L c 
50 for L= Ley. 


Accordingly, any ZL, such as Ly) or Ly, must be obtained by estimating a characteristic of 
the assumed distribution. For reasons described in appendix B, the distribution characterizing 
ball-bearing fatigue life was taken to be the Weibull distribution. In brief, this distribution 
ean be derived by assuming a “weakest-link”’ concept of fatigue strength. In addition, the 
suitability of the Weibull distribution for fatigue life has*been verified in many cases by empirical 
plotting of data 

One method of estimating Ly) or Ls makes use of special probability-plotting paper so 
designed that a theoretical Weibull distribution plots as a straight line, and treats the problem 
as one of straight-line fitting by conventional least squares procedures. However, the procedure 
usually followed does not take into full account the number of bearings that remain intact when 
tests are incomplete, nor the interdependence of successive points. Because of these and other 
limitations, it was decided to use an alternative approach in the estimation of Ly) and Ls for 
each test group (see appendix B). 

To this end, a method was developed that takes into account explicitly the number of 
bearings remaining intact at the termination of a test, and that also possesses several other 
advantages. This method makes use of certain specially determined linear functions of the 
observed failure times (in logarithms), 7;, arranged in order of size. These functions have the 
general form 

T=) €,2; (3) 


-—— 
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As the method makes intimate use of the ordered arrangement of the data, it is termed an 
“order statistics’”’ method. 

The coefficients ¢,; in eq (3) allow great flexibilitv. They have been determined in such 
a manner that the method will have certain desirable objective characteristics, e. ¢., freedom 
from systematic error and a minimum standard error. 


2.2. Evaluation of the Parameters in the Life Formula 


Once the estimates for Ly) and Ly are obtained, it is possible to evaluate the exponent p 
in the life formula. However, in order to make the most efficient use of the given data, it is 
necessary also to estimate the other parameters, f,, a@,, and do. 

The methods for estimating the values of Ly) and Ly) for each test group actually vield 
results for In Ly and In Ly. Thus, taking logarithms * of the life equation (2a) gives 


In L=(p In f.)+ (pa,) In Z+ (pa) In D,—p In P. (4) 
This equation can be written more simply as 


Y =b,-4 bia, + bors i D) 
where 


Y=In L (for either Ly» or Lz), 
bo=p |n f., b;=pa,, b,=paz, b, P, (6) 
2,=In Z, Ze=in D,, %=in P. 


The quantities 7,, 2, and x; depend on the characteristics of the bearing type and test 
conditions, and can be regarded as known exactly. On the other hand, the variable Y, which 
depends on the outcome of the bearing tests, is subject to considerable dispersion. Thus, 
estimates can be found for the parameters bo, b,, b:, and bs, using standard least squares methods 
based on minimizing the sums of squared deviations in the y direction. These methods are 
discussed in detail in appendix C. 

After the parameters ho, b,, b., and 6, are estimated, values for do, a), ad, and p can be 
found from the relations 


f bh, b, 

a nN. bs a “a 
i 

(ly —; p — h. 


It is clear that the values for do, a,, andl a, depend on the value of p. 

The estimates for p and the a’s are subject to some uncertainties because they are based 
on test results, which themselves are subject to considerable variability. Hence with every 
value of Pp and of the a’s calculated from the life data, there is given also an interval of un- 
certainty to indicate its precision. These intervals are “‘95-percent confidence limits.”” * 

A large interval of uncertainty associated with an estimate indicates poor precision; a 
small interval of uncertainty is evidence of high precision. These intervals of uncertainty 
not only reflect the inherent variability of the test data, but are also affected by (a) how well 
the life equation (2a) is the proper functional form for bearing life, and (b) the suitability of 
the data (including the number of test groups) for estimating the parameters in the life formula. 

Further technical details concerning the evaluation of the parameters in the life formula 
are given in appendix C. 

Natural logarithms to the base e are used throughout 


‘ Briefly, confidence intervals describe the compatibility of the observations with an unknown parameter estimated from them; 95-percent 


confidence limits are limits such that on the average, in repeated applications of the same procedure, 95 percent of intervals so calculated will 
contain the unknown true value of the parameter Che confidence limits associated with p are symmetric However, the confidence limits asso 
ciated with the a’s are asymmetric because of the dependence of the a’s on p 
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3. Summary of Analyses 


3.1. Evaluation of Parameter p 


The statistical analysis based on all deep-groove ball-bearing data from companies A, B, 
and © ° yielded the final values for p shown in table 2. The separate values for each of the 
three companies are given in table 3. The intervals of uncertainty specified by + quantities 
in those tables refer to intervals within which, with reasonable assurance, the true value of 
the parameter is located. The fact that all of the intervals of uncertainty exhibit considerable 
overlap shows that the data are consistent with the supposition that all three companies have 
a common value of p for deep-groove bearings. The fact that all the intervals include 3 
indicates that all the estimates of p are consistent with the practice of taking p=3. More- 
over, the value of p for Ly was not significantly different from that for L5p. 


TasBLe 2. Final over-all values of p for deep-groove 
bearings 


p 2.8740. 35 2. 80 +0. 31 


TABLE 3. Individual estimates of p for deep-groove bearings by company 


Number 
Company of test Lio Lso 
groups 


\ 50 3. 00 + 0. 64 3. 05 + 0. 60 
B 148 2.75+ .48 2. 62+ 0. 40 
C l: 3. 12 8S 2. 88 + 1. 02 


The values given for p are based on analyses of all deep-groove ball-bearing data, irrespec- 
tive of bearing type. Hence, the parameter estimates represent “omnibus” values. In order 
to investigate the dependence of the exponent p on bearing type, the data from company B, 
which was made up of three bearing types, were analyzed separately. The results for the 


> 


exponent p are shown in table 4. These results are all compatible with the value p=3. 


TABLE 4. Value of p by bearing type for company B 


Number Value of p 
Type of test 
groups 
/ Ls 
B-1 37 3. 36+ 0. 6S 3. 2340. 47 
B-2 4 2. 65+ 0. 91 2. 13+0. 79 
B-3 17 1. 89+ 1. 28 2. 82+ 1. 10 
Total 148 
The data furnished by company D were too few te be included in the analysis 
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3.2. Evaluation of Parameters f,, a;, and a» 


The computations that give estimates for the exponent p also vield estimates for the quan- 
tities In f., a;, and a. From the relations (6) it is clear that the values for these parameters 
depend on the value for p. Thus, associated with every value of p will be corresponding 
values for In f., a,, and a,. Table 5 summarizes these parameter estimates associated with 
the final values of p. The estimates for ag=In f,, rather than f,, are given here, because this is 
the parameter that arises naturally in the life formula (ef. eq 4). 

The analyses conducted separately for each company resulted in other values than those 
in the previous paragraph for do, a,, a. These results are summarized in table 6. They show 
excellent agreement with the results in table 5, even though the values for p are somewhat 
different. 


TABLE 5 Final values of ao. a;. ao for Ly and L 


Company p Ay Interval of ay Interval of a Interval of 
uncertainty uncertainty uncertainty 
Ly 
2. 87 9 02 4. i, 10. 79 0. 380 0. 454, 1. 201) l.é2 1 51 1. 92 
RB 2. 87 8. 55 7. 98 9. 14 670 0.418. 0. 920) 81 1. 70, 1. 92 
C 2. 87 9. 56 6. 85. 12. 42) . 174 1. 750, 1. 352) 1. 37 0. 09 2. 67 
Lso 
| 
\ 2 80 10. 36 8. S81, 11. OS) 0.015 0. 741. 0. 751) 1. 69 1. 50, 1&8 
B 2 80 9 O05 8. 54, 9. 60) . 695 170, 920 1 91 1. SI 2. 01 
( 2. 80 9 O5 6. 61, 11. 58) 75 . 921. 1. 847 1. 76 0. 60, 2 493 
TABLE 6. Values of ao, a,, a2 for Ly and Ley, based on inde pendent analyses for each company 
Company p a Interval of a Interval of a Interval of 
uncertainty uncertainty uncertainty 
L 
3. 00 8. O7 7. 18, 10. 90 0. 390 0. 507, 1. 249 i. 43 1. 50 1. O4 
B 2. 75 8. 59 7. 99, Y, 24 666 398. 0. 928 1 80 1. 67 1, 92 
C 3. 12 9, 21 4. 29, 11. S4 O41 1. 326, 992 1. 36 0. 49 2. 30 
Lis 
3. 05 10. 13 & 48. 12. 00 0. O72 0. 768 0. 855 ‘#4 1. 50 1 oO} 
B 2. 62 9 15 8. 61, 9. 76 690 156, 22 1. 90 170 2 00 
C 2. 88 8. 93 6. 58, 12. 39 510 1. O55. 1 810 1. 75 0. 66. 3.05 


Similarly, the values for ao, a,, and a), arising from separate analyses made on the three 
types of bearings from company B, resulted in still other estimates for these parameters. Table 
7 summarizes these estimates. These estimates are less precise than the corresponding omnibus 
values given for company B in table 6. This is a consequence of the fact that within a bearing 
tvpe, the quantities Z and J), hardly vary at all. This condition makes the data unsuitable 
for estimating the associated unknown parameters, do, a;, and ay. 
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TABLE 7. Values of ao, a), a2 for Ly and Ly by bearing type (company B) 


Type ay Interval of ay Interval of a2 Interval of 
uncertainty uncertainty uncertainty 
Ly 
B-1 7.25 | ( 5& 21, 939) | 1.07 0.35, 1.76) | 1. 68 (1.27, 2. 08) 
B-2 7. 34 5.54, 9.33) | 1.21 | | . 43, 2.00) | 1. 69 (1. 38, 1. 93) 
B-3 2. 50 7. 95, 16.04) | 3. 70 2.25, 9%. 06) 1. 27 (0. 30, 1. 65) 
L; 
B-1 7. 39 5.92, 8. 92) 1. 23 0. 73, 1. 73) 1. 79 (1. 50, 2. 08) 
B-2 9.00 | ( 7.10, 11.67) | 0. 87 . 05, 1.68) | 1.77 (1. 46, 2.03) 
B-3 1. 03 1.34, 6.53) | 4. 50 ( 1.97, 7.19) | 1.48 (1.21, 1.70) 


3.3. Redetermination of the Estimates for f,. 


The uncertainty intervals associated with estimates for the parameter do=In f, are quite 
large. This is primarily because the uncertainty associated with the estimate of a also de- 
pends on how well the other parameters, a), a), and p, are estimated. Another way to evaluate 
a, Which may result in smaller intervals of uncertainty, is to assume a priori values for a), do, 
and p, and then determine the estimate for a. This procedure was followed by using the 
widely accepted values for the parameters given in [15], namely, a;=2/3, @=1.8, p=3. 

However, if on such a calculation the values assumed for the parameters a,, d,, and p are 
not compatible with the given data, then values of a (or f,) so calculated will not be correct 
determinations for these data. Accordingly, an analysis was made to determine whether the 
parameter values in [15] were compatible with the given data. 

This analvsis showed that these parameter values are compatible with the data, with re- 
spect to all individual companies for rating life Lj, but not for median life Ly. (Company 
A was the only company for which the parameter values are suitable for median Jife.) A fur- 
ther analysis, by bearing tvpe for company B, showed that the above parameter values are 
not suitable for the rating life Ly with respect to B-3-type bearings. 

In the light of this last analvsis, redetermined values of do, taking a,;=2/3, a2=1.8, and 
p—3, are only strictly valid with respect to company A, company B (B-1, B-2), and company 
(' for rating life Ly. These values are summarized in table 8. For convenience, these new ° 
estimates are given for f,=In7'do CXPp do. 


TABLE 8. Values for f, assuming a,;= 2/3, a2=1.8, p=3.0 for Lio 
Company Number of i. Interval of uncer- 
test groups tainty 
\ 50 4,538 | (4, 273, 4, 817) 
B (over-all valu: 148 4, 925 (4, 750, 5, 105) 
BOI 37 1,709 (4,403, 5, 034) 
B-2 O4 5,033 | (4,885, 5, 187) 
B-3* 17 - 
C 12 3, 294 | (3,029, 3, 583) 
D 3 1, 639 (3, 478, 6, 187) 


*Assumed values of parameters a), a2, and p not compatible with test results for bearings of this series, 
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4. Appendix A. Summary of Original Data 


This appendix summarizes in tabular form the worksheets submitted by the American 
Standards Association Subcommittee to the National Bureau of Standards for statistical 
analysis. Separate tables are presented for deep-groove data from companies A, B, C, and D. 
These four tables (A—1 to A-4) are followed by table A—5, which gives a svnopsis of the number 
of test groups and the number of bearings for each company. 

Tables A-1 to A-4 give the size of test group, the values for quantities ?, Z, )),, and the 
estimates ° for Ly, Ly, and the “Weibull slope” e. All of these variables are directly observed 
or specified quantities except for the estimates Lyo, Ls, and e. These last three quantities are 
based on statistical calculations that made use of the results of individual endurance tests. 
These calculations are explained in appendix B. 

The original data, as submitted, contained a few cases where companies tested bearings 
manufactured by other companies. Such test groups are not included in the summary tables, 
as these results confound differences in testing with differences in manufacturing. Therefore 
these test results were not used in any of the analyses. Thus, table A—3, for company C 
omits 4 tests performed on other manufacturers’ bearings; table A—4, for company D, omits 3 
tests. 

The five tables described above are followed by a specimen worksheet * with identifving 
information removed. A sample of Weibull-function coordinate paper is also included. This 
coordinate paper had been used for graphing the results of all the individual endurance tests 
and these graphs had accompanied the worksheets submitted to the Statistical Engineering 


Labora tory. 


Ihe estimates for Ly and Ly are given in milli ition r all companies except company D r fe estimates for D are . 
hours, the same units in which the origina! endurance 
Bearings marked “Omitted”’ were completely e! ted from consideration, as company representatives explained t 
fatigue failures and should not be regarded as part of t test group As aresult, t test [ he case tl pecimen 1 4 
consist of 23 bearings rather than the original number o his type of situation appeared rather infrequently, however 
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TABLE A-1 





Summary ball-bearing data for company A, with computed values for Ly, Ly, and Weibull slope e 

















Record Year Number Z Da e 
No. of in test] Load Number [Ball diam. L Leo Weibull 
test group of balls = slope 
lb in. 
1-1 1936 2h 4240 8 11/16 19. 84.5 1.27 
le 2 1937 20 420 8 11/16 26.2 The2 1.81 
l- 3 1937 1h 420 8 11/16 11,1 68.1 1.0) 
1- 4 1937 19 420 8 11/16 11.8 66.8 1.09 
le 5 1937 18 420 8 11/16 13.5 79 elt 1.06 
le 6 1938 21 2530 9 1/2 5.80 2567 1.27 
le 7 1938 28 420 8 11/16 18.3 bhe7 2.10 
le 8 1938 27 42ho 8 11/16 5.62 7302 e73 
1-9 190 20 420 8 11/16 15.8 82.7 1.Lh 
1-10 190 22 420 8 11/16 8.70 41.6 1.20 
1-11 190 19 42h0 8 11/16 11.6 160 072 
1-12 190 15 19,0 9 7/16 20.6 Tle 1.52 
1-13 | 190 15 190 9 7/16 14.5 88.2 1.04 
1-1) 19,0 15 2536 9 1/2 12.1 33-1 1,87 
1-15 | 190 14 2536 9 1/2 15.1 46.4 1.67 
1-16 | 190 iS 2536 9 1/2 14.0 4326 1.66 
1-18 190 26 420 8 11/16 46.2 110 2.17 
1-19 | 190 14 420 8 11/16 30.0 88.2 1.74 
1-20 | 192 20 420 8 11/16 21.1 S7 el 1.89 
1-21 192 20 420 8 11/16 17.3 4567 1.94 
1-22 192 37 42h0 8 11/16 3725 118 1.64 
1-23 192 36 2h0 8 11/16 20.3 771 1.41 
1.2) 192 32 42h0 8 11/16 4.03 42.5 280 
1-25 191 28 25h 8 17/32 8.38 84.7 281 
1-26 1943 23 3975 8 19/32 1.79 1365 093 
1-27 192 30 4,00 10 5/8 11.7 “Sel 1.39 
1-28 192 31 6920 8 7/8 4.15 15.8 1.41 
1-29 | 1943 30 990 9 5/16 7023 41.0 1.09 
1-30 1943 30 1509 7 7/16 2269 110 1.20 
1-31 | 1943 30 932 7 11/32 9.54 31.6 1.57 
1-32 19h); 26 3180 8 19/32 6.28 2320 1.45 
1-33 | 19kh 29 3180 8 19/32 4.81 21.2 1.27 
1-3) 194 33 860 10 7/8 4.17 12.8 1.68 
1-35 19k) 26 | 14080 8 1-1/4 5 l2 31.6 1.07 
1-36 | 1951 28 190 9 7/16 747 49.5 1.00 
1-37 | 1951 3h 2330 9 7/16 4.80 21.3 1.26 
1-38 | 1951 27 1550 9 7/16 14.8 78.4 1.13 
1-39 | 1951 29 1165 9 7/16 84.9 4,60 1,11 
1-40 | 1951 27 2910 9 7/16 3240 16.5 1.19 
1! 1951 27 3880 9 7/16 1.2) 3.23 1.97 
1-42 1951 26 776 9 7/16 2h1 951 1.37 
1-3 1951 30 | 19750 8 1-3/h 3.01 12.6 1.31 
Lely) 30 2112 8 11/16 89.1 4.86 1,11 
14S . 30 22h 8 11/16 15.2 104 098 
© 
146 a 30 8448 8 11/16 220k, 10.2 1.17 
1-47 30 2112 8 5/8 51.0 376 09k 
1-48 3 30 | 22k 8 5/8 5 026 58.8 ~78 
1-49 30 84,8 8 5/8 883 he 1.09 
1-50 | 19h) 30 422k 8 11/16 14.8 57h 1.39 
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Tasie A-2. Summary ball-bearing data for company B, with computed values for Ly, Ly, and Weibull slope « 

Record Year Number Z D, e 

No. of in test Load Number {Ball diam. yo Ley Weibull 

test group of balls slope 
1b in, 

2-1 19h0 19 570 10 3/16 6.68 13.4 2.72 | 
2= 2 19k 20 570 9 5/16 29.8 70.0 2022 

2= 3 196 23 580 9 1/4 16.3 55 el 1.55 

2= i 19,6 23 580 9 1/4 28.5 69.2 2013 

2= 5 1947 23 580 9 1/4 16.4 49.3 1.71 

2= 6 1943 10 665 9 1/k 10.3 40.1 1.50 

2= 7 1944 10 665 9 1/4 2567 h6.4 3.19 
2= 8 1942 19 580 10 1/4 9.55 3966 1,32 

2= 9 1946 33 620 10 1/4 1729 62.1 1.51 

2-10 19447 15 620 10 1/k 19.9 7302 1.45 

2-11 1947 31 620 10 1/4 12.9 50.4 1.39 

212 194k 19 625 10 1A 19.3 46.2 2.18 | 
2=13 19,1 17 720 10 1/4 11.1 233 2.54 

21); 1946 60 980 11 9/32 15.7 43.5 1.85 | 
2-15 1947 32 980 11 9/32 11.2 38.1 1.54 

2=16 1950 49 600 11 5/16 417 809 2.85 | 
2=17 19h49 60 600 11 5/16 216 709 1.58 

2-18 1943 20 900 11 5/16 35.26 100 1.82 

2=19 1946 67 1220 11 5/16 12.0 42.2 1.50 | 
2=20 19h7 34 1220 11 5/16 8.53 46.6 1.11 

2—=21 190 20 1370 11 5/16 6.77 18.9 1.85 | 
2=22 1950 60 1415 11 5/16 13.5 46.5 1.53 

2=23 1950 60 22)3 11 5/16 2032 8.06 1.51 

2-2) 192 20 720 12 S/16 |} 367 141 1.40 

2=25 1946 55 1300 12 5/16 | 19.0 5702 1.71 

2=27 194) 20 1650 1h 11/32 |} 17.0 | Thek 1.37 

2=28 1946 59 1760 1h 11/32 20.9 5307 2.00 

2=29 1947 34 1760 14 11/32 9256 40.7 1.30 

2=30 1940 20 2010 13 13/32 5049 | 33.3 1.05 

2=31 190 9 2010 13 13/32 1.39 4h.O 054 

2=32 19k) 19 210 13 13/32 9,80 82.7 088 

2=33 1943 i 2630 15 13/32 5219 54.9 280 

2=34 19k2 12 5900 1h 19/32 6.36 17.5 1.86 

2—35 1947 19 5900 1h 19/32 3.68 22.1 1.05 

2=36 1947 20 8070 1, 23/32 83h 2346 1.61 

2=37 192 12 8075 1h 23/32 6.78 36.4 1.12 

2-39 190 2 720 8 9/32 18.2 56.9 1.66 

2-10 190 2h 720 8 9/32 22.8 56.2 2.09 

2-1 1940 25 720 8 9/32 3.99 15.6 1.38 

2—))3 190 25 720 8 9/32 el! 28.5 1.36 

Qalihs 190 25 720 8 9/32 12.5 26.) 2051 

2-5 19h1 25 72 8 9/3 18.8 48.7 1.98 

2-16 1941 23 720 8 9/32 21.5 5302 2.08 

2=)9 1943 20 900 9 5/16 30.1 92.3 1.66 

2650 | 19kh | 18 | 00 9 ial isco | h7.6 1.63 
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Tape A-2. Summary ball-bearing data for company B, with computed values for Ly, Lo, and Weibull slope e—Continued 

















Record |} Year Number yA D, a 
No. of in test Load Number Ball diam. No Ley Weibull 
test group of balls slope 
j lb in. 
2-51 195 27 94,0 9 5/16 17.5 52.8 1.71 
2=52 1947 34 94,0 9 5/16 Wek 656 1.24 
2-53 1938 10 1180 9 5/16 8.76 22-1 2.04 
2-5), 195 30 1580 9 3/8 12.1 43.3 1.47 
2-55 1947 33 1580 9 3/8 17.2 6406 1.42 
2=56 1948 8 1580 9 3/8 10.7 3426 1.41 
| 2=57 19h5 31 2160 9 7/16 10.9 3726 1.52 
2-58 197 30 2160 9 7/16 12.7 5307 1.30 
2-59 1938 9 2200 9 7/16 3.73 4305 e77 
2=60 1947 30 2,80 9 7/16 16.6 78.3 1.21 
2-61 1950 ho 130 9 15/32 180 275 elals 
2-62 | 1937 19 1660 10 7/16 852 23 1.86 
2-63 191 19 1700 9 15/32 S7el 230 1.35 
| 2—6), 1939 2h 21,80 9 15/32 15.7 55.8 1.48 
2=65 1939 25 2,80 9 15/32 2761 97.8 1.47 
| 2=66 1939 23 2,80 3 15/32 21.7 122 1.09 
2-67 1939 28 2,80 9 15/32 13.2 42.3 1.62 
2-68 1939 28 2,80 2 15/32 35.8 145 1.35 
| 2-69 | 1939 20 2,80 9 15/32 12.7 3407 1.87 
2-70 19h) 20 2,80 9 15/32 10,1 2728 1.87 
| 2=71 1945 20 2,80 9 15/32 8.83 34.3 1.39 
2=72 1938 10 2,80 9 15/32 16.5 6003 1.45 
2=73 192 11 2,80 9 15/32 1729 65.8 1.45 
2-7) 1943 10 2,80 9 15/32 15.67 63el 1.35 
2=75 1943 20 2,80 9 15/32 10.8 42.1 1.38 
2-76 | 194k 18 21,80 9 15/32 14.2 3909 1.83 
2-77 | 19bh 18 280 9 15/32 19.0 67.8 1.48 
2=78 19)); 18 2,80 9 15/32 16.3 577 1.49 
2=79 194 20 2,80 9 15/32 2.93 18.0 1.04 
2=80 19 20 21,80 9 15/32 5.69 2564 1.26 
2-81 19h 28 2,80 9 15/32 9.54 3909 1.32 
2-82 1944 22 280 9 15/32 12.6 5567 1.27 
2-83 194); 23 21,80 9 15/32 510 37-5 094 
2-8 19h 18 2,80 9 15/32 16.0 5307 1.56 
2-85 194) 20 2,80 9 15/32 1.98 22.1 78 
2=86 1945 20 21,80 9 15/32 5.65 28.8 1.16 
2-87 195 20 2,80 9 15/32 12.8 4306 1.58 
2-88 1945 20 2,80 9 15/32 9.8) 3223 1.59 
2-89 | 19h5 20 214,80 9 15/32 12.1 43.0 1.48 
2-90 195 20 2,80 9 15/32 548 40.8 29k 
2-91 195 20 2,80 9 15/32 6.64 2563 1 
-— 19he 32 2,80 9 15/32 1329 41.9 1.70 
293 19,6 35 280 9 15/32 9,02 Look | 
2-94 194.6 34 2,80 9 15/32 11.0 49-2 1.2 
0-95 19h:7 31 21,80 9 15/32 14.5 7306 1.16 
2-96 19h 9 2,80 5 15/32 5.91 372 1,02 
9 19h) 10 2180 3 15/32 18.1 40.5 2.33 
2-98 1945 10 2,80 9 15/32 17.1 533 1.65 
2-99 | 19h5 10 280 9 15/32 32.6 61.8 2095 
2=100 | 1945 10 2,80 9 15/32 24.1 66.2 1.87 
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TaBLe A-2. Summary ball-bearing data for company B, with computed values for Ly, Lso, and Weibull slope e—Continued 












































Record Year Number Z De 3 
No. of in test | Load | Number Ball diam. Lio L Weibull 
test group of balls 50 slope 
1b in. 

2-101 1945 20 2,80 9 15/32 36.1 71.6 2075 
2-102 1946 20 2,80 9 15/32 63.3 104 3.82 
2-103 19,6 12 2,80 9 15/32 paren 59.0 1.33 
2-10) 1946 11 2,80 9 15/32 15el 92.9 1.04 
2-105 1945 10 2,80 9 15/32 18.8 390k 2255 
2-106 1950 12 2,80 9 15/32 5.63 37 1.0) 
2-107 1950 12 24,80 9 15/32 7.23 34.5 1.21 
2-108 1951 30 21,80 9 15/32 16.7 71.8 1.29 
2-109 1951 63 2,80 9 15/32 26.5 90.3 1.54 
2-110 1950 23 24,80 9 15/32 8.35 49.1 1.06 
2-111 1943 19 3250 9 15/32 3-79 9-30 2.10 
2-112 1937 10 34,70 10 7/16 9205 36.6 1.35 
2-113 194 20 4,000 9 15/32 2.98 7035 2.08 
2-11) 1943 19 2300 10 15/32 22.5 730k 1.59 
2<115 1938 10 2730 10 15/32 3.82 31.7 289 
2-116 194.6 22 2660 10 17/32 6.55 20.8 1.63 
2-117 19h 2 2250 11 15/32 1725 6403 1.45 
2-118 19443 16 2300 11 15/32 61.7 152 2.10 
2-119 1945 48 28,0 11 15/32 18.6 4207 2.27 
2-120 1947 28 28,0 11 15/32 21.6 66.3 1.68 
2-121 1947 8 280 11 15/32 11.9 39.1 1.59 
2-122 1948 8 2840 1 15/32 1309 50.6 1.6 
2-123 1943 19 3200 1 15/32 7280 3301 1.30 
2-12) 1944 28 4,000 11 15/32 3055 139 1.38 
2=125 1943 19 4,000 1 15/32 90 23h 2.06 
2126 1947 23 6350 11 11/16 4.76 22.7 1,21 
2-127 194) 20 12000 11 1-1/16 3223 9,86 1.69 
2-128 19h 20 12000 11 1-1/16 2.62 9.52 1.446 
2-129 1944 9 12700 8 1-1/2 7.89 3927 1.17 
2-130 199 18 16500 nu 1-1/16 4.93 20.4 1.33 
2-131 1950 20 16500 11 1-1/16 6.26 16.2 1.98 
2-132 1938 8 565 7 5/16 3703 103 1.85 
2-133 1944 20 900 7 5/16 14.0 38.6 1.86 
2-13) 1938 10 1650 8 13/32 30.3 87.6 1.77 
2-135 19k 20 2250 8 15/32 2567 71.2 1.85 
2-136 1943 20 2300 8 15/32 10.5 6.k 1.07 
2=137 19k) 19 3200 8 15/32 10.3 2h.1 2.21 
2-138 19k) 19 000 8 15/32 4.56 12.9 1.81 
2-139 1937 10 1710 8 17/32 251 27h 079 
2-140 1938 9 2360 8 17/32 48.8 264 1.12 
2-11 1937 10 2680 8 17/32 7253 60.7 290 
2-142 1937 i 3850 8 17/32 14.9 62.6 1.32 
2-13 1947 21 7760 8 29/32 4.57 430k 8, 
2—Lhh 1943 12 9550 8 1-1/16 3.90 40.7 -80 
2-145 1947 21 9750 8 1-1/16 15.5 79 04 1.16 
2=1)6 1948 16 11,00 8 1-3/16 10.2 43-9 1.29 
2-147 1948 20 114,00 8 1-3/16 4.71 16.9 1.48 
2-18 1947 18 11420 8 1-3/16 10.1 34.2 1.55 
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TABLE A-3. Summary ball-bearing data for company C, with computed values for Li, Lo, and Weibull slope e 
Record Year Number Z De € 
No, of in test | Load Number | Ball diam, Lio Leo Weibull 
test group of balls slope 
lb in, 
3-1 192 9; 1580 7 9/16 16.9 64.8 1.40 
3 2 199 29 79 7 9/16 211 729 1.52 
3= 3 1949 35 1185 7 9/16 Thelk 287 1.40 
3- 4 19)¢ 29 1600 9 1/2 9.62 40.1 1.32 
3- 5 195 kK 1600 6 15/32 11.9 66.3 1.10 
3- 6 1943 9 2275 7 17/32 13.8 58.0 1.31 
3-7 196 13 2510 . 15/32 2.38 11.3 1,21 
” ~ 1946 12 250 O 15/32 2238 11.5 1.19 
3e § 19h 12 1580 7 9/16 8.75 62.2 0.96 
3-10 1949 12 1580 7 9/16 2507 113 1.27 
3-11 197 2h 1600 9 1/2 1h.5 113 0.92 
3=12 19,9 12 610 g 5/16 26.8 6546 2.10 
TaspLeE A-4, Summary ball-bearing data for company D, with computed values for Ly, Ls, and Weibull slope e 
— 
Recorc Year Number Z Da e 
No, of in test Load Jumber | Ball dian, dl Le Weibull 
test group of balls slope 
lb in, 
u-1 1946 19 1750 7 7/16 159 963 1.05 
le3 | 19s. | 56 1750 9 7/26 113 582 | 1615 



































* life estimates are in hours. 
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TasBLe A-—5. Summary of test groups of ball-bearing data 


Number of test Total number of 
Company groups bearings in test 
group 
A 50 1, 259 
B , 148 3, 289 
Type B-1 37 ‘ 
Type B-2 94 . 
Type B-3 17 aa 
C 12 291 
D* 3 109 
Total (all companies) 213 4, 948 





*These data were not used in the main analyses, 
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SPECIMEN WORKSHEET 



















































































Reference No. Table Ordered According to 
Bearing Mfg. by Endurance Life 
Bearing Tested by 
Date of Test 2-26-16 
Bearing No. Brg. Endurance Type of 
Load 5cC RL. No. Mill.Revs. Failure Rema rks 
Speed 2CCO_r.p.m. 1 7 20 e 
Lubrication: Type__ Jet Oil = ee eo 
Frequency___ e a0 0 all 
Ball No. and Dia. 9 = 4" 19 11.52 I.R. 
Contact Angle 0 9 eo Rell 
Groove Radius: Inner Ring _ 51.6% ll eae val] 
Outer Ring____§3.0% 15 hehe Ball 
Number of Rows 1 12 51.¢h ball 
Bore 20 mm. 20 51.96 Ball 
0.D. 42 mm. 18 ch el I.R. 
Lot Size 25 Taken on 23 13 55.56 I.R. 
1 67.8C Bal] 
et 
Bearing temperature measured on outer a ae x an Parra 
ring at point of maxjmum load h 68.6, Bell ee 
Material: Type 6 68.54, L.Bor 
Source 5 Le ao ies 
Rockwell Hardness of: - 3) 1D > any 
Inner Ring 63,5 17 93.12 Ball 
Outer Ring flos 7 08.64, I.R. 
Balls 23. -105.12 I.R. 
24 106.84 —®> Disc. 
oe 
21 127.92 Ball 
Ball Failure 13 52% 3 128 4, O.R. 
Inner Ring Failure 5 20% lL 172 LO => Dise 
Outer Ring Failure} LZ ae eat i 





Test life in 10° revolutions: 

















Median 6s. 
Mean v2 * 
B-10 29. 
Slope of Curve 2.23 
Test No. 3183 
Lot 7] 








—_—— _ th ccc 
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SAMPLE OF WEIBULL FUNCTION COORDINATE PAPER 








BEARINGS TESTED- PERCENT 
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5. Appendix B. Evaluation of Ly», L», and Weibull Slope e, by Using Order 
Statistics for Censored Data 


This is a technical appendix that gives the mathematical basis for estimating, for each 
test group, the values of Ly) and Ls for use in the regression analysis discussed in appendix C, 
and also the Weibull slope « 


5.1. Weibull Distribution 
a. Characteristics 


As noted in the text, the basic assumption for estimating Ly, Lo, and e for each test group 
was that the probability distribution of fatigue lives of individual bearings could be represented 
by a “Weibull distribution.” * This means that the observed fatigue lives (number of revolu- 
tions) of all the bearings in a test group of, say, n bearings constitute a random sample of n 
independent observations from a distribution whose cumulative (from above) distribution 
function (hereafter denoted by edf) is ° 


S(L)= Prob {life>L} 
exp |[—(L/a)‘}, 0<L<o, (Bl) 


where a and ¢ are the two parameters to be fitted. They are related to Lijp and Ls by eq (B2a) 
below. The function S(L) is also termed the “survivorship” function. This distribution is 
one of three limiting types to which the distribution of the smallest member of a sample, under 
general conditions, tends as the sample size is increased indefinitely. (Another type is dis- 
cussed in the following section.) This matter was first studied chiefly by Fisher and Tippett 
[5], and for this reason the type (B1) is sometimes referred to as Fisher-Tippett type III for 
smallest values. 

There are both theoretical and practical reasons for choosing the Weibull distribution 
(B1) as the underlying probability distribution for fatigue life. 

Theoretical. Here it is assumed that fatigue is an ‘“‘extreme-value”’ phenomenon, related 
in some manner to the strength at the weakest point in the material under stress. The theo- 
retical reasoning that proceeds from this assumption is mentioned by a number of authors, 
and is given explicitly, for example, by Freudenthal and Gumbel in [6, p. 316 to 318]. It 
leads precisely to the form (B1) (see eq (2.9) in [6]). It is recognized that this statistical 
assumption has not received universal acceptance. This paper is, however, not concerned 
with the relative merits of various statistical theories of fatigue, but merely with consequences 
of a reasonable choice from among them. 

Practical. Application of the Weibull distribution received extensive attention by W. 
Weibull in [19], where he showed that a distribution of the general type (B1) represented certain 
fatigue-life data quite satisfactorily. In addition, inspection of the special ‘‘Weibull’’ plots 
accompanying the worksheets suggests that many can be fitted satisfactorily by a straight 
line representing a Weibull distribution, as explained below. 

The manner in which these graphs are constructed is described by Weibull in [19]. A 
sample of Weibull-function coordinate paper used for this purpose is included in appendix A. 
The essence of the method is that eq (B1) may be converted, by taking logarithms twice, into 


e(In L)—(e In a)=In{[In(1/S)), (B2) 


where “In” denotes the natural logarithm (base denoted by e) and S=S(L). From eq (B1) 
and the definitions of Ly and Ly, when L=Ly», S(L)=.90; and when L=Ly, S(L)=.50. 
These values substituted in eq (B2) give 


* So named for W. Weibull (cf. [18], p. 16 ff.), who is considered to be one of the first to study it extensively 
* The use of a continuous instead of discrete probability distribution will introduce no appreciable error. 
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e(In Ly) —e(In a) =In{In(1/.90) |= —2.25037 
(B2a) 
e(In Ls) —e(n a) =In{In(1/.50)|= —0.36651, 


the values on the right-hand side being obtained from [17, table 2]. These are the relatioa- 
ships between the parameters a, ¢, Ly, and Lz. The right-hand numerical values will later 
be denoted by 9, Y¥.9, respectively. Equation (B2) may be written 


ex—a’=y, 


where 
r=In L, a’=e In a, y=In{In(1/S)]. (B3) 


The variables x, ¥ correspond to the two scales shown on the Weibull-function coordinate paper 
in appendix A. The variable z, with unrestricted values, corresponds to the horizontal scale 
“Bearing life,” having a logarithmic seale. The variable y is represented through the per- 
centage surviving, S, or rather through the (vertical) scale for “bearings tested—percent’”’ = 
percent failed '°=1—S=P, which can vary only between 0 and 1. This seale also has non- 
uniform graduations, given by the iterated logarithm in (B3). 

The Weibull distribution is thus seen to be equivalent to a straight-line relationship, with 
“Weibull slope” ¢, between the logarithm of fatigue life and an associated quantity y depend- 
ing only on its relative rank when the fatigue lives are arranged in ascending order. Thus, 
goodness of fit of the straight line (B3) is equivalent to goodness of fit of a Weibull distribu- 
tion to the fatigue lives Z of an individual test group. In fact, one common method of sta- 
tistical analysis of fatigue-life data (Freudenthal and Gumbel [6]) depends upon the use of 
the classical method of least squares for fitting this straight line. This method is, however, 
subject to certain limitations described below. Instead, an alternative method, presented in 
the following sections, is preferred that fits the distribution of z=In Z directly by use of order 
statistics. 


b. Limitations of Fitting by Least Squares 


In the classical method of least squares for fitting the straight-line relationship (B3) to 
a test group of ball-bearing data, pairs of values (2, y,, i=1, . .., m, are required. The 
values of z=In L are obtained from the given data. However, the variable y, measured 
through the percentage failing, P=1—S, presents difficulties. The problem of how to plot 
P is known as the problem of “plotting position.”’ 

It seems clear that the values, ?;, of the plotted variable, ?, must somehow be related to 
the rank order of the bearings as they fail. A natural choice is the percentage failing: P=f/n, 
where f is the rank order of failure in a test group of n. This is not advisable for reasons dis- 
cussed at length by Gumbel in [9, p. 14], where he advocates the plotting position f/(n+1)." 
Other workers take different positions, and the question of plotting position must be regarded 
as still unsettled. 

A second difficulty with the use of least squares is that as usually used it fails to take 
adequate account of the number of items remaining intact (“runouts’”’) in the incompleted 
tests. Asa final point, it is to be noted that the successive plotted points are not independent, 
as they represent the observed lives in increasing order. A correct use of least squares pro- 
cedures would have to take into account all the intercorrelations, which is not done in the 
usual application of the “method of least squares.’”” The method of order statistics described 
in section 3 has the advantage of avoiding the above limitations of the least squares method. 


lhe symbol P as used here should not be confused with the same symbol for load used in the life formula. In any event, the meaning 


will be clear from the context 
' This plotting position was also used by Weibull in [18] (ef. eq (72) and the vertical scale in figures 3 and 4 therein 


289 








5.2. The Extreme-Value Distribution 
a. Relation to Weibull Distribution 


The preceding section indicates that logarithms of lives, rather than lives themselves, are 
the natural units in which to carry out the analysis. This idea has also been adopted by those 
who do not use the Weibull distribution, either because they are unaware of its existence or 
because they do not feel it fits their data. 

If the Weibull distribution is adopted for fatigue life, 1, then the variate, r=In Z, has 


the nonnormal cumulative distribution function 


G(r)= Prob! In (life) >z Prob /life>¢ 


S(é)=exp(—a“e”) exp| ge- man") ocr ow, 
This may be written 
G(r) =(y) =exp(—e), (B5) 
where 
y= (r—u)/B, oC Z< ©. (B6) 
and 
u=Ina, B=I1/e (B7) 


are the two parameters. The distribution, &(y), considered as a distribution of the “reduced 
variable,” y, has standardized parameters u=0, 8=1, and is called the ‘reduced distribution.” 

The form (B5) is another of the three asymptotic distributions of extreme values, some- 
times designated as Fisher-Tippett type I for smallest values. This distribution has been 
studied extensively, chiefly by E. J. Gumbel (e. g., [7, 8, 9]). In this paper, the term ‘“‘extreme- 
value distr.bution” will be given to the distribution of smallest values (unless otherwise speci- 
fied), although this name is frequently given to the largest-values case. 

From the above discussion, it is apparent that methods pertinent to the type I extreme- 
value distribution (B5) are appropriate. For this purpose there is available a mathematical 
approach recently developed by one of the authors of this report, and described in detail in [13]. 


b. Characteristics 


A description of the extreme-value distribution (B5), together with an interpretation of its 
parameters in terms of life estimates (or rather their logarithms), is essential to an understand- 
ing of the application of the method of order statistics in this paper. It will be seen that the 
problem of estimating life is equivalent to that of estimating the parameters u and 8. 

The parameters of the extreme-value distribution (B5) are depicted in figure 1 (page 291). 
The quantity u is the position of the mode or highest point of the (frequency) distribution. 
The quantity 8 is a scale parameter, analogous to the standard deviation, ¢, in the case of the 
normal distribution. In fact, 8 is \6/m (about %) times the standard deviation of the extreme- 
value distribution. 

Although the two parameters, u, 8, completely specify the distribution, it is very useful 
to introduce related quantities of the form 


t=u-+ By, (BS) 


which are linear combinations of parameters u and 8 and may thus also be regarded as param- 
eters when known values are later assigned to y. Introduction of ¢ makes it possible to esti- 
mate uw and 8 simultaneously. Thus if ¢ can be obtained as a+ by with v and 6 known and y 
arbitrary, then the values «=a, B=) can be read off at once. 

The parameter ¢ has another highly important meaning. In figure 1 the area F under the 
distribution to the right of the ordinate erected at ¢ represents the probability that a value 





2 Cf. Freudenthal and Gumbel [6], eq (2.5), (2.6), (2.8), (2.9). 


290 
































yp =— In(—InF) 
W/) 





k-BY- 








x 
© 
i // 
: Y, 
S X19 =INLig=U+BY_ 96 
a. 
o 
rT Y 90> ~ 2.25037 
za 
F =.90 
> Y, YY) 
: WY 
> 4 
> t—— BY 59 ——4 
WW te= In Lig U 
oO 
X50=!NL5o=U+BY 50 yy Y; 
Y.50 = ~ 0-3665! F=.50 
te=InL so t BY.s50 











FIGurReE 1, 


x=In L 


General form of extreme-value distribution (for smallest values) showing relationship of 
parameters tr, Bo In Lio, and a——In Ly, to u and B. 
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larger than ¢ will occur. Thus ¢ is a function of F and may be written fy, as shown; it is desig- 
nated the “upper 100F-percentage point” of the distribution. For example, if /=.90, then 
tt o represents a value of r=In L, which will be exceeded by 90 percent of the population. 
This is associated with rating life 1,, (life exceeded by 90 percent of bearings) by the relation 


$co=Zi9= In Le, (B9) 
where z represents life in logarithmic units. Similarly, for median life, 

to=Za=In Lo. (B10) 
Since the ¢’s are regarded as parameters of the distribution, so also are x;) and 259, and therefore 
Ly and Ly. These are not, of course, all independent. 


In general, we have the percentage point fy, which, expressed in terms of the original param- 
ta] t z 
eters u and 8, may be written in the form (B8): 


tp=u+ Byp, (BS8a) 

where y is a quantity depending only on the probability 7, determined as follows. We have 
from (BS8a) 

Yr — (te —u)/B, (B11) 


i. @., Yr is the value of (x—1u)/8 when « takes the value fp. But by definition of the probability F, 
in view of (B5), (B6), and (B11), 


F=Prob{r>t, G (te) =P(yp) = exp(—e*). (B12) 


Thus, solving for ye, we obtain 


Yr—\|n(—In F). (B13) 


This is the reduced variable corresponding to the probability F’, and may be obtained by a simple 
change in sign from table 2 of [17], which tabulates the function 


In(—In #,), 


where ®,, a probability, takes on values from 0 to 1. Thus, 


for F=.90, yp=— 2.25037; B 
(B14 
for F=.50, yy 0.36651. 


y? 


The above discussion shows that both s;) and 25) (rating and median lives in logarithmic 
units) may be determined once the general percentage point (BS8a) is estimated by giving the 
two particular values (B14) to yp. 


c. Conversion From Largest to Smallest Values 


The methods and numerical results developed in [13] were for problems, such as maximum 
gust-loads on airplanes, that required the distribution of /argest sample values. In order to 
adapt this material to the distribution of smallest values (B5) required here, the relationships 
of symmetry involved in the reversal of direction must be examined with care. To avoid 
confusion, it is necessary to use subscripts Z and S to distinguish between quantities related 
to the largest-value distribution from those related to the smallest-values case. No generality 
is lost by use of reduced variates. Thus, in (B5), # will be replaced by the reduced variate 
y, and, for simplicity, the symbol G(y) will be used instead of &(y): 


G(y) = ®(y) =exp(—e’). (B15) 
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From this, the (“‘cumulative from above’’) distribution of smallest values is 
Prob{ Ys>y} =Gs(y) =exp(— 2), o<my<oa, (B16) 


where }’s denotes the reduced smallest value. The corresponding distribution of largest 


values is (see Gumbel [9, eq (1), p. 21)) 
Prob{ Y,>y H,(y)=1—Prob{ Y¥.<y} =1—exp(—e") = 1—G,(—y), —7<y<oa, (B17) 
from (B16). 
The corresponding relation for the density functions is obtained by differentiation, with 
gs(y) = Gs(y), and hy (y)=Hz(y): 
gs(y) =hi(—y). (B18) 


Hence the two distributions are merely mirror images of each other. The moments of the 


distributions are related as follows: 


ns = Ely") | y‘gslydy | (—y')*hily’)dy=(—1)' nec. (B19) 
Thus, the means differ in sign and the variances are identical: 
ns——Y Vins (B20) 
» FF. » 
=——s1- (B21) 
) 


These values are given, for example, in [9, p. 23, eq (3.27)]. 

Finally, we need the relationships between moments of the order statistics for the two 
distributions. As the smallest-value distribution is a reversal of the largest-value distribution, 
it is natural to reverse the arrangement of the order statistics as well. This gives simpler 
results. Thus we are interested in the ith order statistic in the series 

(S): yiz>yz2>... DY>.-- Vay (B22 
where the parent distribution is that of smallest values. Primes will be used as a reminder 
that the order is descending, not ascending. Thus in tables B-2 and B-3 the absence of primes 
indicates that the order statistics are in increasing order. 

(B22) is the analogue of the series 


(iL: mime... SRS .« + Se (B23) 


of order statistics for the largest-value parent distribution. Whenever a distinction is neces- 
sary the subscripts S or Z will be used with the y’s. 

From (B18) it may seem intuitively (and may be justified rigorously) that the distribu- 
tions and moments of the order statistics follow the same symmetry relationships as the parent 
distributions, namely, 

E.(y,)=(—1)"E, (y?) 
Es(yiy)= Eviyy,) 
> (B24) 


Os j =0'(Y) OLY) =9F, 3 





Os. 69 =FS(Yi, Ys) =O (Ys, Y)) =x 513+ J 


In other words, the even moments remain the same; the odd moments change only in sign. 
The above development shows that the numerical results for moments of order statistics 
previously obtained in [13] for the largest-value case can be used here for smallest values 


without any substantive change. 
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5.3. Method of Order Statistics for Censored Samples 
a. For Small Samples 


Consider an independent random sample of x items from the distribution of smallest 
values, of which only the & smallest values can be observed. In view of the preceding dis- 
cussion, it is desirable in the theoretical development to deal with the order statistics in descend- 
ing order: 

(A ST A 


> > Ree? are 2 (B25) 


where the parentheses denote the (n—s) (largest) unobservable values, and the remaining k 
values are known. ‘This arrangement materially simplifies the exposition. Primes will again 
be used to denote descending order to distinguish from ascending order, which will occur in the 
later parts of this section. 

From the & known values it is desired to determine an estimator 


Te. x=WiFn—eai tpt Wodta-teot «1. tWete, KK, (B26) 
(i. e., the weights w) of the general parameter, 
tp—u-+ Byp, (B27) 


of the extreme-value population (B5), such that 7” in (B26) is (1) unbiased and (2) of minimum 
variance. Mathematically, this means that 


ECT’) =f,, (B28) 
where / denotes mathematical expectation, and 


Var (7’)=a minimum, (B29) 


subject to the above condition. 
From (B6), 
r=u-+ By, (B30) 


where y is the reduced variable and x« the observed variable. From this the following relations 
for the order statistics 2; and y, are apparent: 


2; u+By;, j=n—k+1, n—k+2,.. .,n, (B31) 
x, >] ae ya 7 ae (B32) 
Vo-2412Yn-2422>- - +2 Ve (B33) 

E5(2',) =u+ BEs(y}). (B34) 


The values s(y) may be obtained with the aid of the table in [14]. This table gives the 
values of /,(y)) where the order statistics, y), are in descending order (as indicated by the 
prime). The means needed in (B34) are obtained from (B24) and the evident symmetry 
relations 

Ely?) 1)" Fp (yn 141) 
as 
Es(yi) Ex (Yn—s+1)- (B35) 


Reference [14] gives the values of /,(y)) for r=1(1)min(n,26), n=1(1)10(5)60(10) 100, 
From (B28) and (B31), 


E(T"\=>) w;, [u+BE(y),-24))=tr=u+ Bye. (B36) 
j=1 


k 
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This is required to be an identity for all values of the parameters u, 8. Equating their co- 
efficients gives the two conditions on the weights, u 





k 
=| , 


Wi=1, DLE (yn—n+3)) W3 


Yr, 
j=1 


S 


k 
2 (B37) 


where the numerical values E(y;,_;,;) may be obtained from [14] as already indicated. 
For the variance condition (B29), we have, in view of (B26), 


k k 
o(T")= 2 w7?o(tp—v4 Jt 2 DW Wyo(Ty—e41) Cn—e4s) (B38) 
j=! j=l] i=] 
i-X) 
Var(7")=[2w)?on 24 jFEZ'WiWion_e 41, nee 1P=VE" BP 


(B39) 


minimum subject to (B28). 
Use of Lagrange multipliers in the same manner as in [13, pp. 50-52] gives, after differentiation, 
the conditions on the weights: 


; k 

,2 ; . ’ .y e J ’ 
On—k4 Ww; -S On—k+ii,n—k4 W; T r T ulvy, k4 0,7 m4 ee ey k. 
i=! 

ixj 


, 


(B40) 


For each fixed value of kK<n there are k linear equations which, with the two in (B37), form a 
simultaneous system of (k+2) equations in the (k+2) unknowns, wj, wo, . . ., wi, 4, wu. ~The 
values of \ and uw are useful as a check, because, if (B40) is multiplied by w, and summed for 
j=1,2,. . ., k, the result is, in view of (B37) and (B39), 


V et d T uy p (). 
The minimum value, V{, 


(B41) 
k, min» W ill be denoted by Qy. ». 


In general, there will be a set of (k+2) linear equations to solve for each k=2, 
(1) Casek=n. Fork 


ee 
n, the matrix of coefficients and right-hand “constant terms’”’ of 
(B40) and (B37) is the (n+2) by (n+3) matrix 








xa ae oe Tin l Ey; o* 
o> O59 ae Ton l Ey; 0 
A° : F (B42) 
wae o,9 ate. Son l Ey, 0 
] l l 0 0 ] 
Ky, Ey Ey, OO 0 Ye |. 
The ordinary (n- 


2) by (n+2) matrix of coefficients, without the constant terms, will be de- 
noted by A,. If T, denotes the vector column of constant terms, then 


) 
AS=[A,/T,], 


(B43) 
and the linear system of (n+2) equations may be denoted 
A,W,=T,, (B44) 
where W% denotes the column vector of the (x 


») , , 
+2) unknowns W), W», . 


i, 
295 








The coefficients of the unknowns in (B44) involve the means F(y}), already discussed, and 
covariances o,,. These values are given in table B-1 for n=2 to 6. The oj, were computed 
by the method developed in [12]. Table B-1 also indicates how the moments for the largest 
values case can be obtained simply from those shown. 

The (n+ 2) solutions of (B44) are all expressible linearly in terms of the components of T,. 
Thus the solutions all take the form 


w=a’+ bye, j=1,2, .. 8 
A=ce,+diy, (B45) 
u=e,+dy, 


Substituting these w) in (B39) gives an expression of the form 


(, V iw = A+ 2 Baye t+ Cry? (B46) 


dn. min 
The quantities a}, 6; for the weights w), for n=2 to 6, are shown in table B-2. The coefficients 
A, B,, C, of Qh.,, and the values of Q},, evaluated at F=.90, .50, corresponding to Lyo, Ls. 
respectively, are given in table B-3. 

Calculations were limited to n=6 in this paper, in view of the diminishing returns in 
“efficiency” (see below) for increasingly larger amounts of computing. Methods suitable for 
larger values of n are discussed later. 

Table B-3 shows that as sample size increases from n=2 to 6 (in the case k=n), the 
variance diminishes for the percentage-point parameters ft for F=.90 and .50, 1. e., fe=sy—In 
Ly and te=F5——=In Ly. This is a common characteristic of the behavior of estimators for 
increasing sample size. Another method whereby estimators may be compared is through 
their efficiency. 

Efficiency is a measure intended to provide a convenient standard of comparison for 
estimators. This is done for two estimators to be compared by dividing the variance of each 
into a theoretical “smallest” variance, Q,,, known as the ““Cramér-Rao lower bound.” Further 
details in the case of complete samples where k=n, as here, may be found in [13, p. 14 and 
15); values of Ure are also indicated in this reference in table IIT (a). 

Table B-4 shows the efficiency values so obtained, for the case k=n, n=2 to 6, as regards 
the order-statistics estimators for the parameters sj In Lio, Zeo=In Leo. 

These values show that for x, the efficiency, starting with under 70 percent for n=2, 
increases rapidly until 89 percent, out of a possible maximum 100 percent, is reached for n=6. 
A 90-percent-efficient estimator is generally considered to be quite good. As regards £5, 
the efficiency is well above 95 percent for all the values of n, and for n=6 exceeds 99 percent. 
In view of results of this nature, and because of the increasingly heavy computations necessary, 
calculations were not carried beyond n=6 in [13]. 

The above applies to estimation of the parameters sj, and zo, which it will be recalled are 
the logarithms of the actual life estimates Ly, Ls. It is believed that efficiency of the method 
of order statistics in obtaining estimates of actual life Ly), Ls is probably reasonable, in view 
of its high efficiency in estimating the logarithms, Y10, Tso. 

(2) Case k<n. For the case k<n, the procedure is very similar. One starts with a 
(k+2) by (k+3) order matrix A? derived from A’ in (B42) by striking out the first (n—hk) 
rows and columns. One proceeds in this manner for k=n—1, n—2, ete., until when k=2 
the matrix becomes 


[ 1,n—1 oh ln l Ey, I 07] 
oe Ton ] ky}, 0 
A} (B47) 
I l 0 0 l 
| Ey’ , ky; 0 0 ved, 








representing a set of 4 equations in 4 unknowns. 
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The resulting weights w; and variances (;,,, were obtained in similar manner to those 
for k=n in (B45) and (B46). These, it will be recalled, are primed quantities, associated with 
descending order of the order statistics. Because the observations, xz, for successive failures 
naturally occur in ascending order, it is more useful for actual application, in contrast to theo- 
retical development, to tabulate the weights and covariances for the order statistics in ascending 
order. This has been done in table B-2, giving the weights w,;—a;+byp, and in table B-3, 
giving the variances (,,».—A+2By,+ Cy; for the estimators 7T,,, formed with the above 
weights. These variances are also evaluated for the parameters z;>=In Ly, %=In Ly. The 
relationships of these unprimed quantities to the primed ones of the previous theoretical 
development is merely a reversal of the order throughout, as indicated by subscripts: i. e., every 
a’, is changed to the corresponding a,_;,; and similarly for 6; and wi. The variances Q may 
be shown to remain unchanged. 

b. Extension to Larger Samples 


Samples of more than six items are broken up into independent samples of 6 with a re- 
mainder subgroup, if necessary, of from 2 to5items. Because the endurance data were arranged 
in increasing order of life, independent random subgroups could not be obtained by simply 
taking groups of 6 in the (numerical) order in which they appeared on the worksheets. It was 
therefore first necessary to randomize the endurance lives on each data worksheet. This was 
accomplished by use of random numbers that were generated in the electronic computer (the 
SEAC) as needed. 

Such artificial randomization is not desirable when it can be avoided, because the results 
of the caleulations are then not unique, but may depend to a limited degree on the particular 
set of random numbers used." It is therefore recommended that when the bearings in a test 
group are to be simultaneously run on a battery of fatigue-testing machines, the individual 
bearings should be recorded in advance in some more or less natural order independent of the 
order in which failure takes place in the course of the test. Natural order might be order of 
manufacture, order of testing, ete. 

In the present investigation, each subgroup was treated as a random sample by the methods 
already developed for size 6 or less. That is, a “subestimator’’ was calculated for each sub- 
group and the results averaged to produce an over-all sample estimator. 

An estimator, both for the individual subgroup and for the over-all sample, was obtained 


for each of the four population quantities: 
u, B, to=Zo=In Lyo=Uty.8, tso=So=In Lyo=Ut ¥ Af. (B48) 


For subgroups, these four parameter estimates are given by (caret denotes “estimate of’’) 
| | a 


A I ‘ A a m A 
u= > a,2;, B=>, 5.2; Tyo= Ty, C10)=U+Yy. 0B, 50 T, (50)=t-4 Y.50B, (B49) 
l i=] 
where ¥ » 2.25037, Ys 0.3665l,andz,SmS .. . Sa,2SkSnS6, are the logarithms 


of the actual observed lives in a subgroup arranged in ascending order, and the a; and 6, are 
read directly from table B-2. For the over-all sample estimator, the subestimators 7, , are 
merely averaged. 

For later use (appendix C) the variance of the over-all estimator, 7, and its relation to 
sample size will be considered here. Consider first the case of a complete sample, where no 
intact bearings are present because the test is run to completion. Let n be the sample size; 


then there are two cases, according as (1) n<6, or (2) n>6. 


Chis effect can be reduced somewhat by making a duplicate run and averaging the results, as was done in this study. 
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Figure 2. Relationship of variances Qy and Qs) to sample 
size n for n 2 to 6 (logarithmic scale in each direction). 


Ow is variance of estimator of z In Zao 
Os is variance of estimator of 2 In Lso 
All 9’s are in units of 8? 


(1) n<6. Table B-3 gives the numerical variances, Yio and Qs, for n=2 to 6. These 
values are plotted in figure 2 on double-logarithmic paper. The values for (so (right-hand 
scale) are seen to lie on a straight line of slope negative unity. This shows that at least in this 
case, variance is inversely proportional to sample size. For the other case, (jo, a straight line 
also gives a reasonably good fit, and its slope appears to differ only a little from —1. Hence 
the underlined statement is approximately applicable here too. 

(2) n>6. Ifa sample of size n>6 is broken into equal subgroups (of size 6, for example) 
there will generally be a remainder of size less than 6. The preceding development, when 
suitably modified, shows that for large n the influence of this remainder is small compared to 
the remaining bulk of the sample and thus the rule in question holds approximately in this case. 
Agreement with the rule is less close for a few cases of moderate n, but for simplicity the inverse 
relationship will be taken as a reasonable rule of thumb in all cases for the over-all purposes 
of analysis. 

Two complete runs were made on the SEAC for each of the 213 test groups of data, and 
the two results were averaged for each group, giving values of the averages 


u, B, - In Lio, te In Les. (B50) 


From these, the values of Ly and Lo were obtained from a table of exponentials and the 
Weibull slope e=1/8 obtained as a consequence of formula (B5). An example showing the 
steps in calculation of Lyo, Lo, and e is discussed below. 

Because this investigation represents probably the most extensive mass fitting of the 
Weibull distribution made to date, a tabulation of the 213 values of the parameter e will be of 
considerable value to future applications of this distribution. This is shown in table B-6. 
The corresponding histogram is given in figure 3. Particular items of interest are 


mean of ¢ (all 4 companies)= 1.51, 
median ¢ (all 4 companies) = 1.43. 


Note that 50 percent of the values are in the interval e=1.17 to 1.74. 
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iuRE 3. Histogram for the 213 estimates of parameter e for companies A, B, C, and D combined. 


_ 


Data in table Bi 


c. Worked Example 


The example that will be given to illustrate the foregoing procedures will be the one that 
was worked out as a “‘test problem” for the SEAC before using the full set of data. The test 
group of bearings selected for this purpose was No. 1—1 in table A-1, for company A. The test 
group consisted of 24 ball bearings, of which 4 remained intact when the test was discontinued. 
The details of the computation for obtaining values of Ly, Ls0, and Weibull slope e from the 
test group of data are contained in table B-5 and described in the steps below. 

The endurance lives in observed (increasing) order are listed in column (1). The arrows 
indicate the four “run-outs,” or “‘intacts,’’ whose testing was discontinued at the number of 
million revolutions indicated. All that is known about these four bearings is that their fatigue 
lives exceeded the values shown. 

Step 1. Randomization. The order of endurance lives in column (1) was randomized by 
use of a set of random numbers generated in the SEAC as part of the computation work. The 
result is shown in column (2) of table B—5. 

Step 2. Subgroups. The lives in randomized order were divided, as shown by the lines of 
separation, into subgroups of size n=6, the maximum size for which the order-statistics weights 
had been computed."* Each subgroup was then prepared for the application of the order- 
statistics method by rearranging in increasing order (column (3)). Natural logarithms were 
then taken as in column (4). 

Step 3. Weights. Each subgroup was regarded as consisting of k actual observations out 
of a censored sample of n. It happened here that n was 6 for every subgroup; k took the values 
6,5, 6,3. These values are shown in the subscripts of 7, , written in column (3), and they 
determined the weights a; and 5, to be selected from table B-2. These weights are represented 
in columns (5) and (6). 

Step 4. Cross-products. The cross-products 


k 


ah =) = 
(> py ms >, O42; 
i=] i=1 


were then evaluated and placed as shown for each of the subgroups. 


‘ As sample size 24 is an exact multiple of n=6, it so happened that there was no ‘‘remainder subgroup” in this case. This will not usually 
be true, but the procedure is identical for other values of n, differing merely in the numerical weights to be used. 
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Step 5. Estimates. 


the two columns (5) and (6), and denoted by 7; 
statistics estimates of the two parameters u and 8 of the extreme-value distribution that 
represents the underlying Weibull distribution. 


A . 
The reciprocal of 8 yields the Weibull slope « 
The following logarithmic life estimates were 
a 


a 
8, respectively. 


A OS 
u, T, 


1.32497. 


. - . . 
the estimates u and 8, using the given values of ¥ 9, ¥ 50: 


rating life Lio 
median life Le.» 


These three values (A), (B), 


a 
L109 In Ine u 


In Lep=t 


antilog 


a 
(Zio) (base e) 


. - 
antilog (Teo) (base €) 


~ 


ys 2.982305 


A 
25037 B 
= - 
0.36651 B=—4.404120 
19.2333 million rev 


81.7872 million rev. 





A simple arithmetic average of the four values was taken for each of 


These are the order- 


(A) 


given by following linear combinations of 


and (CC) represent the outcome of the calculation. 


In the full-scale computing program, calculations were carried out by the SEAC to a larger 


number of places than is shown in the table for presentation purposes. 
the number of places shown here should be adequate. 


The values Ly, Leo, and 


In gene 


ral, however, 
( shown here 


differ slightly from those recorded in table A-1 because the latter represent averages of two 


separate runs. 


TaBLe B-1. Means, variances, and covariances of order statistics y; in samples of n from the reduced extreme- 
value distribution Gly exp ey), 2lo6 
For distribution of largest values, y;<yo< . <n 
For distribution of smallest values, y/> y’ A 
Means* Variances and covariances,* ¢/.=¢..—¢e a 
” 7 Es y. ) 
E, Yi 
} l =Z } 5 ) } f= & 
au J 0. 11598 152.0. 68402 804 0. 48045 301 
- i @ 1. 27036 285 1. 64493 407 
| 0. 40361 359 14849 7960. 30137 1440. 24375 S10 
5 2 . 45943 263 65852 235 . 54629 438 
2 1. 67582 795 1. 64493 407 
l 0.57351 263 . 34402 417 . 22455 3440. 17903 4540. 15388 918 
{ 2 LOGOS 352 $1553 113) . 33720 966 20271 ISS 
3 . 81278 175 . 65180 236 . 57432 356 
t 1. 96351 003 1. 64493 407 
l 0. 69016 715 . 28486 447 IS202 536 14358 737 0. 12257 8650. LO9OL 329 
2 10689 454 . 80849 748 . 24676 731 . 21226 644 IS967 383 
5 3 $2555 O61 10598 292 . 35267 O72 . 31716 095 
t 1. O7093 582 . 64907 319 . SS99L 519 
5 2. 18665 358 1. 64493 407 
| 0. 77729 368 . 24658 20 15496 74 12121 61 . 10291 64 0. OOL16 19 
2 . 25453 448 . 24854 56 19670 62 16806 28 . 14945 32 
¢ 3 ISS838 534 20761 59 . 25616 60 . 22887 90 
, { . 66271 588 10185 52 . 36145 55 
5 1. 27504 579 . 64769 96 
6 2. 36897 513 


*The means are for smallest values (denoted by subscript S 
; are the same for both. 


order of y’s The o 
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0. OS285 42 
13619 10 
20925 46 
33204 51 
5VOSS OF 

1. 64493 41 


- for largest values, change all signs and reverse 























TABLE B-2. 


distribution (smallest values) from a censored sample of n 


bo 





to 


tbo 


Weights w; for the order-statistics estimator T,,,; for the parameter tr 


a 


a 


b 


a 
lh 


OS36269 


. 7213475 
. 3777001 
. 8221012 


OS79664 


. 3747251 


. 7063194 


S690 149 


. O80 1057 


1143997 


07 13800 


2487965 


. 9598627 
. 8962840 
. 2101141 


1345419 


0153832 


. 2730342 


0583502 
1844826 
1655650 
9141358 


. 3153968 


1466018 


. 0865378 
. 2858647 
. 0057311 


2015451 


. DASS6OBY 


1458072 


. 916373 


. 7213475 


. 3777001 
. 8221012 


. 2557135 
. 2558160 


. T0638 194 
. 8690149 


. 06043 16 


0. 
. 2034315 
. 3886492 
. 0280534 
. 2654739 
. 0465729 


. 3258576 


1536799 
2239192 


. 9598627 
. 8962840 
. 086023 1 
. 3642463 
. 0519642 
. 2499429 


LO88236 
1816564 


9141558 


1972753 


. 083522 


221 
1495332 


- +Wele,Wj=a;4+ biUp 


<r 


rr, k=2ton 


0. 6563201 


6305411 


. 0196741 
. 7402573 


. 2639426 
. 0859035 


. 2961372 
. TY85882 


1520750 
1491094 
1676091 
1304534 


. 5188283 


83525 10 


. 0649390 


301 


IS5S8756 
1002523 
1536040 
1210527 


1267277 


0. 


rs 


5109975 


. 5586192 


. 8115440 
. 6720865 


. 2462831 
. 0065354 


. 0496521 
. 7372142 


1722784 


. 0645894 


1656192 


. 0731937 


u+Byp of the extreme-value 
2 to 6, where only the k smallest values are known 


0. 4189341 
. 5031278 


. 6751653 
. 6170118 
. 03859868 


Ts 


0. 3554481 
1592751 





2 


TasBLe B-3. Variance Q,.,8? of order-statistics estimator T,.,., given in table B-2, and its numerical values 
Qn.x (10) =Qio, Qn.x(0) =Qso for estimators of parameters t.»=2o=—lnLy, t.5.=—2s=—lnLy, respectively, for a 
censored sample of n=2 to 6 

Variances in units of 8? 


ask 1+2Byp+Cy,2 


< k 2ton 
n k 1 B ( ( Ose 
2 2 0. 6595467 0. 0643216 0. 7118574 3. 975015 0. 708021 
3 (2 . 9160886 1682465 8183654 2. 952920 . 682735 
; | 3 1028637 . 0247719 . 3447117 2. 260033 . 467327 
‘2 1. 3340189 7720208 . 8670220 2. 250056 _ 884572 
| 3 0. 4331573 1180273 . 3922328 1. 888278 . 399329 
} . 2934587 . 0346903 . 2252828 1. 590460 349150 
2 1. 7891720 1. 0115594 . 8950462 1. 769068 1. 167910 
. | 3. 0. 5293953 0. 2353740 . 4168155 1. 580861 0. 412852 
- 1 2918142 . 0385708 . 2537913 1. 403458 _ 297633 
5 . 2313953 . 03839905 . 1666472 1. 228307 . 278697 
(2) 2. 2440055 1. 2082248 . 9132926 1. 431164 lL. 481035 
3 0.6529409 0. 3332488 _ 4321160 1. 341381 0. 466709 
6 | . 3237185 1020223 . 2697162 1. 230430 . 285165 
5 _ 2236063 . 0105329 . 1861069 1. 118677 240885 
6 . 1911738 0313731 . 1319601 1. OOO644 231897 
*For O..y, 2.25037: for Yep 0.36651, 


Tas_e B-4. Efficiency of order-statistics es- 
timator of logarithmic life xy=Iln Ly and 
t~=—ln L for complete samples (k=n) of 


size 2 too 


i: fficiency (in percent 
with re spect to 


J In L Ix In L 


2 67. 2 97. 4 
3 78.8 O98. 3 
! S4. 0 98. 7 
5 87. 0 Os. 9 
6 So. 0 QQ. | 
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TABLE B-—5. 


data 


Test group 1-1, company A 


Endurance (million revolutions 


Natural loga- 


Observed Randomized Ascending rithms x; 
order order order within 
subgroups 
6. 0 119. 0 27. § x,=3. 31419 
8. 6 138. 0 69. 0 2=4. 23411 
17.8 146. 0 119. 0 z3—=4. 77912 
18. 0 151. 0 138. 0 X% 4. 92725 
27.5 21.5 146. 0 r5;— 4. 98361 
33. 5 69. 0 151. 0 x6= 5. 01728 
T 3. 
0. 5 150. 0)— 6. 0 x,=1. 79176 
1. 5 & 6 8. 6 z 2. 15176 
69. 0 1.5 1.5 r 3. 94158 
74.0 89. 0 89. 0 ry= 4. 48864 
74. 0 109. 0 109. 0 r 1 61935 
So. 0 6. 0 150. 0 r 
T,,. 
109. 0 74.0 18. 0 r, = 2. 89037 
118. 0 ISL. O 33. 5 I 3. 51155 
119.0 141.0 14.0 I 4. 30407 
138. 0 18. 0 118. 0 ry;=4. 77068 
141. 0 33.5 141. 0 r 1. 94876 
144. 0 144. 0 144. 0 r 1. 96981 
T., 
146. 0 17.8 17.8 r,;= 2. 87920 
150. 0 153. 0) 90. 5 r= 3. 92197 
151. 0 153. 0) 74.0 I 1. 30407 
153. 0 153. 0) 153. 0) ry) 
153. 0 50. 5 153. 0) r:) 
153. 0) 74.0 153. 0)— Xe) 
T, 
Over-all estimator 
SUMMARY 
T,=a=4, 6807388 7,—8=0. 754735 
/ 9 2.25037 Y 50 
Zo=In Li=t+ y 06 = 2.982305. 
to=In Ly = t+ y 508= 4.404120, 
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Law, 


k 
Laz, 
l 


k 
Faz, 
l 


k 
Laz, 
l 


f 


0.36651 


Example showing use of order statistics method of com puting values of Lio, Ls, and e from endurance 


5 6 


Weights 


a b; 

0. 048867 0. 145807 
. 083522 . 149533 
. 121053 . 126728 
. 165619 . 073194 
. 225491 . 035987 
. 355448 . 459275 
4. 817310 Sb 0. 400992 

0. 005731 | 0. 201543 
. 046573 . 197275 
. 100252 . 153604 
. 172278 . 064589 


. 675165 . 617012 


k 
Lh.zi 


1 


. 213655 


1. 446363 


Weights are same as 
for first subgroup 


_ 


. 62808 1 


619440 
446602 
. 388649 
. 835251 


Eb.2, 0). 
0. 315397 0. 
. 203432 
1. 518828 


_ 


1. S31L197 


_ 


784853 


k 
Lbiz; ( 3 
1 


“I 


. 680738 0. 754735 


B= 1. 32497 


Ly= 19.7333 


Ls = 81.7872 








Tarte B-6. Tabulation of estimated values of Weibull slopes, e, for the 213 test groups of companies A, B, C, and D 


Company 


Weibull slope, ¢ Total 
\ B Cc D 
0.50 to 0.74 2 l 3 
.75 to .99 6 1] 2 l 20 
1.00 to 1.24 16 22 3 2 13 
1.25 to 1.49 11 10 5) 56 
1.50 to 1.7 7 32 l 10 
1.75 to 1.99 6 17 23 
2.00 to 2.24 2 12 l 15 
2.25 to 2.49 2 2 
2.50 to 2.74 } 1 
2.75 to 2.99 } } 
3.00 to 3.24 l | 
3.25 to 3.49 
3.50 to 3.94 
3.75 to 3.99 l l 
1.00 to 4.26 
1.25 to 4.49 | l 
Total 50 148 12 3 213 

Mean 1. 33 1. 60 1. 31 1. O4 1. 51 
Median 1. 27 1. 49 1. 30 1. O06 1. 43 
50% interval* 1. O7 to 1. 59 1. 27 to 1. 82 1. O8 to 1. 45 0. 94 to 1. 16 1.17 to 1. 74 


*First quartile to third quartile. 


6. Appendix C. Evaluation and Analysis of the Unknown Parameters in the 
Life Equation with Respect to Companies and Bearing Types 


6.1. Summary 


I B= P 
ss P b 


expresses the dependence of fatigue-life ZL on the design characteristics of the deep-groove 
bearing (Z, D,), the bearing load P, and the “workmanship factor’ f,. This appendix out- 
lines the statistical methods that were used (a) to determine “best’’ empirical values for the 
parameters f,, a, @2, and p of this life formula, (b) to derive the associated intervals of uncer- 


. 


Equation (2a) of the main text, 


tainty, and (c) to answer various questions about the values of these parameters, from the 
basic endurance data furnished by the ASA Subcommittee, which are summarized in appen- 
dix A. These methods are applied separately in each case to the rating life Ly and median 
life Lso values derived from the endurance test data as described in appendix B. 

Section 6.2 summarizes the application of the statistical methods used to determine best 
empirical values and intervals of uncertainty for the parameters f,, a), d:, and p in the case of 
deep-groove bearings. Sections 6.3 to 6.6 outline the statistical analyses employed to answer 
various questions about the values of f,, a, a2, and p for the ball bearings of companies A, B, 
and C. In particular, section 6.3 gives the analysis employed to determine whether values 
of these four parameters are the same for the bearings of the three companies. This analysis 
is carried out separately for rating life Lj) and median life Ly, and the postulated ‘complete 
between-companies homogeneity” is not supported in either instance. Section 6.4 gives the 
analysis appropriate to determining whether the data are consistent with the supposition 
that the value of p is the same for the three companies (regardless how the values of the other 
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parameters may differ); this analysis is applied to the Zy>) and Lso data, with an affirmative 
conclusion in both instances. It is concluded further that the data are consistent with the 
supposition that this common value for p is equal to 3. 

Section 6.5 presents the analysis employed to determine whether the values of the param- 
eters f., @, @,, and p are the same for the three types of deep-groove bearings (B-1, B-2, and 
B-3) for which the data from company B were available. A negative conclusion is reached. 
These data are then reanalyzed to determine whether they are consistent with the supposi- 
tion that p=3 for each of the three types, regardless of differences in the other parameters. 
An affirmative conclusion is reached in this case. 

Section 6.6 concerns (a) the extent to which the Lj) and Ls values are consistent with the 
supposition that the values of @,, a, and p are 2/3, 1.8, and 3, respectively, as given in [15]; 
and (b), the determination of more precise values for f, in those cases in which the foregoing 


supposition is supported. 


6.2. Determination of ‘‘Best’’ Values for the Parameters and Their Associated Intervals of 
Uncertainty 


As shown in section 2.2 of the main text, if natural logarithms are taken of both sides 
of the life equation (eq (2a)), the resulting equation expresses the logarithm of rating life Lyo 
(or median life Ls.) as a linear function of the logarithms of the characteristics of the bearing 
(Z.D,), and the bearing load P, with coefficients that are simple functions of the ‘workmanship 
factor” f, and the exponents, @), @:, and p; 1. e., 


Y=6b)+ 6,2, + boret bars, (C1) 
where 
Y=In L ) 
m—In Z ; 
(C2) 
r=In D, 
zzin P 








and 
bo=p In f, Pao | 
b,=pa, : 
. (C3) 
b, Paez 
b,=—p J 


are unknown constants to be estimated from the data. 

The variables x,, x:, and z; are fixed variates. Their values are uniquely determined by 
the design of the bearing and the bearing loads that are used in the tests. The variable Y, on the 
other hand, denotes the mean values of In Ly, or In Ly, for the population of all bearings with 
characteristics x, and zo, tested at load zs. 

In the practical situation Y is never known, but must be determined from the results 
of endurance tests. The methods used for obtaining such estimates of Y from endurance-test 
data are given in appendix B. To distinguish Y from an empirical estimate of it, the estimate 
will be denoted by the lower-case letter ¥. 

Generally speaking, an estimate y is a random variable, having a probability distribution 
that depends on z;, 2, and z;. We assume that the mean of this distribution is Y= Y (2), 22, 2s), 
and that its dispersion, or more precisely, its variance, is inversely proportional to the rumber 
of bearings w in the test group from which the estimate y=y (2%, 2%, 2; w) is derived (ef. 
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appendix B); that is, 
mean of y=Y 
a (C4) 


variance of y=— 
a 


where o? denotes some positive corstant, and w is the number of bearings in the test group. 

The statistical methods used to estimate the unknown parameters bo, 6,, 2, and 4; from the 
data are termed regression techniques. The books by Anderson and Bancroft [1], Dixon and 
Massey [3], Hald [10], Kempthorne [11], Mood [16], and Wilks [20] give extensive discussion 
of these techniques. For completeness, some of the techniques and rationale of regression 
analysis bearing on the work embodied in this report are summarized below. More detailed 
discussions can be found in the above references. 

Estimation. The problem of estimating the unknown parameters in the life equation can 
be stated as follows: Given independent observations (Ya; Lia; La, L3a} Wa) from n test groups 
(a=1,2,. . .,n) where y, is the estimate of In Ly or In Ly, Xe, 22a are the logarithms of the 
bearing characteristics Z and D),, x5. is the logarithm of the load, and w, is the number of 
individual bearings tested for the ath test group; required to estimate the values of the param- 
eters bo, b;, b., and 63, in eq (C1) using some optimum method of estimation. 

Estimates for the 6; (¢=0,1,2,3) that are free of systematic error and have smaller variances 
than any other linear unbiased estimates are obtained by minimizing the quadratic form, 


n 
@ - > > : Wal Ya by bite boIoa 6,7,.)°, ( ‘D) 


a l 


with respect to each of the 6, (¢=0,1,2,3). 
The resulting normal equations defining the parameter estimates can be written in the 


form, © 
a a A a 7 
oo Do T lo, T Ago» T (log bs Jo 


A “~ a A 
059+ 4416, +-12b2+-A1 35 “Gi . 
> (C6) 


a a : a a 
299+ A215; +22b2+-do3b5 = J2 





A a“ a a 
3059+ 31; +-g2b2+ A335 = 35 J 


where 
n 
a; » ® Wakial ja i); =-0,1,2,3. (C7) 
a=! 
n 
’=>) Wal iala 3==0,1 2.3. (C8) 
a=1 


and Z 41 for all a. 
If one defines the residual sum of squares by 


n 3 A 2 
S=> > W. | ye bare | , 


a=l i=0 


which also can be written in the alternative form 


n 3 A 
S=>5 way2— Do bigv (C9) 
a I 1=0 
then the optimum estimate of o? is 
Ss 
= . (C10) 
n—4 


' A caret (A) is used here to distinguish the fact tnat the solutions of the normal equations are estimated values of the parameters and not 
the parameters themselves 
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The quantity (n—4) is the rank of the quadratic form S, and is termed the degrees of freedom 
associated with S. ing 

The above method of estimation does not depend on ya having a particular assumed 
probability distribution. All that is necessary to specify about the probability distribution of 
Ya is that it possess a finite mean and variance. The solution of the normal equations (C6) 
has the property that the weighted sum of squares of the deviations about y, will be a minimum. 
This is a so-called “‘least-square”’ property of the solutions, but it is only a consequence of the 
method and is not the justification for using this method of estimation. The justification for the 
method is that this is the only one that results in minimum variance unbiased estimates for the 
6, (¢=0, 1, 2, 3). 

Values for 6; (¢=0, 1, 2,3), and hence for a; and p, can be obtained by using either the 
In Ly or In Ly values for ya. In all cases where the a; and p have been obtained for rating 
life Ly, another set of parameters have also been calculated for median life Lo. 

Intervals of uncertainty and inferences. The methods so far discussed for finding estimates 
of unknown parameters need no assumption as to the form of the underlying probability dis- 
tribution of y.. However, something more must be assumed about the distribution of y. if 
(a) one wishes to place an interval about an estimate of a parameter that will include the 
“true” (or population) value of the parameter with given assurance, or (b) if one desires to 
make inferences about the parameters of the life equation for the population from which the 
bearings are a sample. 

Although the endurance lives for individual bearings may follow a Weibull distribution, 
the distribution of y. will not be of this form. However, the estimate y, (cf. appendix B) is an 
average of several independent estimates, each based on linear functions of six or less order 
statistics. Hence, by the central limit theorem, the distribution of the estimate y, will be 
approximated by a normal distribution when n is large (ef. Cramér [2, p. 213]). The statistical 
tests of significance used in this report are not greatly affected by moderate departures from 
normality. Therefore, for making all inferences, it will be further assumed that the estimates 
Ya follow a normal distribution. 

The intervals of uncertainty calculated for each parameter are 95-percent confidence 
limits (equivalent to the usual ‘‘two-sigma”’ limits), which were referred to in the main text. 
Confidence limits for the parameter p= —6; can be calculated by using conventional methods. 
However the confidence limits for ao, a,, and a, are somewhat more complicated. The method 
used here, for this purpose, is sometimes referred to as Fieller’s theorem (cf. Fieller [4]). 

In order to make inferences about the parameters with respect to the different companies 
or bearing types, certain statistical tests of significance were used in this paper. These are all 
based on a test statistic /’, termed the variance ratio or F-ratio, which takes the form 


7 Q, V; 
; Q: vy" 


where Q, and Q, are quadratic forms calculated from the data and », », are the respective ranks 
of the quadratic forms. The explicit expressions for Q,; and Q, depend upon the particular 
hypothesis being tested. The subsequent sections that employ a variance-ratio statistic also 
give the explicit expressions for the two quadratic forms. 

If the hypothesis being tested is true, then the calculated variance ratio will deviate from 
unity in accordance with its tabulated distribution. However, if the hypothesis is false, then 
the variance ratio will be substantially greater than unity, and the “more false” the stated 
hypothesis, the larger the value for the variance ratio. Thus when the hypothesis tested is 
false, this will be detected by an abnormally large F-ratio. In order to objectively determine 
whether a calculated variance ratio is significantly greater than unity, one selects from tables of 
the variance-ratio distribution a critical value of F, such that there is only a small probability 
of the caleulated variance ratio exceeding the critical value from purely chance causes. The 
critical value for F’ used for all variance-ratio tests, in this paper, has been selected so that there 
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is only a probability of .05 of its being exceeded by a calculated variance ratio from purely 
chance causes. This critical value will be denoted by F'9;(%),~). More extensive discussions of 
these procedures are given in Kempthorne [11, chap. 5). 

All statistical tests of significance to test relevant hypotheses have been carried out both 
for rating life Lj and median life Ly. 


6.3. Analysis to Determine Whether Companies Have Common Values for All the Parameters 
in the Life Equation 


This section deals with the details of estimating the values of the parameters in the life 
equation for each company. Furthermore, a statistical analysis is made to determine whether 
the companies have common values for all the parameters in the life equation. The F-ratio 
(eq (C16)), which is used to test this hypothesis, is obtained from the following procedures: 
A single set of parameters, 6, ((=0, 1,2,3), is obtained by first fitting all the data, irrespective 
of company, to the logarithmic life equation, and then calculating the resulting residual sums 
of squares S (eq (C14)) having 206 degrees of freedom. If the hypothesis of common values 
for all the parameters is not true, then a better fit to the data can be made by fitting the life 
equation separately to each company. These calculations result in the individual residual sums 
of squares S,, S:., and S; (eq (C13)) having 46, 144, and 8 degrees of freedom, respectively. 
Thus the total residual sums of squares (S,+S.+-S;) will have 46+ 144+8=198 degrees of 
freedom. Then the difference between S and (S,+S,+-S;), i. e., {S—(S,+S.+S;)}, is also a 
quadratic form having 206—198=8 degrees of freedom. If a substantially better fit was ob- 
tained by fitting a separate life equation to the data for each company, as compared to a single 
life equation, the difference between the two residual sums of squares {S— (S,+S,+S;)} will be 
large. To determine whether this difference is statistically significant the variance ratio 
(C16) is employed. 

Mathematical formulation. It will be convenient to adopt the following notation: Let the 
superscript (w)=1, 2, 3 refer to companies A, B, and C, respectively. Also, for each company, 
let Xe ({=0,1,2,3) refer to the estimates of the parameters in eq (C1), and aj, g}” denote 
the sums of cross products defined in eq (C7) and (C8)"*. Then the normal equations that give 
the parameter estimates for the uth company are 


S* ai b\* qi" 1=0,1,2.3, (C11) 


and the estimates for the parameters a; (¢=0, 1, 2) and p, in the life equation, are obtained from 
the relations given by eq (C3). These results are summarized in tables 6 and 3, respectively, in 
the main text. 

The hypothesis that all parameters in the life equation are the same for each company is 


equivalent to the hypothesis that 
b b\?) —$§ i=0, 1,2,3. (C12) 


Define the residual sum of squares for the uth company by 


Su= 25 WvaYiia— Dy Oi” gi" u=1,2,3 (C13) 
and let 
4 Ne 3 A ; 
S= D5 Do CuaYia— Dy 5191, K=2 9", (C14) 


a=! 1=0 u=1 


where 6; (¢=0,1,2,3) are the estimates obtained from the solutions of the normal equations 


As it is only the ratio of the numbers of observations that is important for weighting, the weights w used in the calculation of the sums of 
cross products in all the analyses have been taken as integral multiples of 5; e. g., if the number of bearings in a test group is 26 (say), then w=5 
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without regard to company differences, i. e., 


[Sai | => 9  i=0,1,2,3. (C15) 
j=0 u=1 u=|] 


Then the variance ratio for testing the hypothesis given by eq (C15) is 


, (S—S,—S.—S;)/8 
k ee he (C16) 
CS, T So T S3), LYS 
and the critical F value is F’»; (8, 198)=1.98. 

Table C-1 summarizes the values for S, (u=1,2,3) and the other quantities needed to 
evaluate eq (C16). 


Paste C-1. Synopsis of calculations to determine whether companies have common values for all 
parameters in life equation 


Lio Lo 
Degrees e 
of 
freedom Mean Mean 
Sy square Ss square 
Companies combined 206 146. 684633 348. 285324 
; 16 102. 552212 89. 734079 
B 144 285. 121962 201. 935797 
C 8 6. 160858 8. 357822 
Sum 198 393. 835032 1. 9891 300. 027698 1. 5153 
Difference 8 52. 849601 6. 6062 18. 257626 6. 0322 
; 6. 6062 6. 0322 
Fr gq (C16 i - 3. 32 : : 3. 98 
nn "=Too1 = ° F=1 5158° 


The values of the F-ratio (eq (C16)) calculated from the quantities in table C-1 are 
Ly: F=3.32 
Ly: F=3.98 


These calculated F values are both larger than the critical value, Fo; (8, 198)=1.98; actually 
the probability is less than .0005 of having an F-ratio as large as those above from purely 
chance causes. ‘Thus from the above statistical tests of significance, the conclusion is reached 
that the three companies do not have common values for all of the parameters in the life equation. 


6.4. Analysis to Determine Whether Companies Have a Common Value for the Exponent p 


The previous analysis resulted in the conclusion that the parameters in the life equation 
are different for each company. However, this does not exclude the possibility that all com- 
panies may have a common p, even though the a; (<=0,1,2) differ from company to company. 
This section discusses the analysis made to determine whether the companies have a common 
value for the exponent p. The analysis given here consisted of the following procedure. First, 
the logarithmic life equation (C17), having a common value of p, but allowing the a; to vary for 
each company, was fitted to all the data, and the resulting residual sum of squares S’ (eq 
(C20)), having 200 degrees of freedom, was calculated. The total residual sum of squares 
from fitting the life equation separately to each company (allowing p to vary in addition to the 
a,) is given by (S,+S,+S;) having 198 degrees of freedom (cf. section 6.3 of this appendix). 
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Then the reduction in the residual sum of squares achieved by using a different exponent p for 


each company is{S’— (S,+S,+-S;)} having 200—198=2 degrees of freedom. To test whether 
this reduction in the residual sums of squares is statistically significant, the variance ratio 
(C21) is employed. 

Mathematical formulation. The logarithmic life equation, having a common value for 
the exponent p, can be written for the ath test group in the uth company as 


Y pe Oe + bi ri + OM re + bg2i” a=1,2,..., %,; #=1, 2, 3. (C17) 


i x - -a s] ja 


Note that although each company has the same parameter 6; in the above equation, the para- 
meters 6}, |”, and 6” are different for each of the three companies. Thus there fare 10 
different parameters, i. e., 6} (¢=0,1,2; u=1,2,3) and d,, to be estimated from the data. 

The normal equations for estimating these parameters are 


Sr ai’b a‘® 0s= q\“ s=xf, 1,2: a=1,3,3 
a = 
-— (C18) 
SS dial } 
a ) a dn q 
" l 0 
where 
; 3 
, x " 
@u= >, as a=), 9° 
u | u=!1 


Thus the set of equations given by (C18) is a system of 10 linear equations in 10 unknowns. 
Once the solutions are obtained, the estimates for a,(7=0,1,2) in the life equation are calculated 
from the relationship 


7. =—- t=O, 1, 2; s=1, 2, 3. (C19) 
b, 


The residual sums of squares (denoted by S’) associated with fitting the life formula (C17) to 
the data is given by 


; n 3 4 
1 = a” 
ta » thay te. ~ . 2a A . ~, 
Ss a 2 WuaYua— 25 2 5\" q bss. (C20) 
u la l u 1 i=0 


Then, to test the hypothesis that the companies have a common value of p, regardless of 
the values for the other parameters in the life equation, the variance ratio 


(S’— S,— S,— S;)/2 


° : : (C21) 
(S,+ So0+S83)/198 
having 2 and 198 degrees of freedom is used. The critical / value is F.9;(2, 198) =3.04. 
The values of S’ for both Ly and Le» are 
Lie. S’=393.272847 d. f.=200 
| * )» )) 
Lo: S’=301.687871 d. f.=200, 


and the caleulated "’ F’ values (using eq (C21)) give 


Le FF 0.141 


i 


Lo: F=0.548. J 


From theory, the calculated value for the F-ratios can never be negative rhe reason for the negative value of F for Lio is that the value 
for the numerator of eq (C21) is only accurate numerically to one decimal place on account of round-off errors arising from the solution of the 
normal equations (C18). Thus, if the hypothesis of a common p value is true, then the F-ratio will nof be large and round-off errors may affect 
the resulting calculation. Alternatively, if the null hypothesis is false, then the calculated F-ratio will be larger then 3.04 and the round-off error 
should be of no consequence. 


310 





— eer ee 




















Because both variance ratios are smaller than the critical value FP»; (2, 198)=3.04, the 
conclusion drawn from this statistical analysis is that the data support the hypothesis of a 
common value of p for the three companies. This holds both for rating life (1,9) and median 
life (Ly). The values for the common p are given in table 2 of the main text. The values 
for the remaining constants in the life equation a,(¢=0,1,2) are found from the relations (C19). 
These results are summarized in table 5 of the main text. 


6.5. Analysis to Determine Whether the Three Bearing Types From Company B Have Con- 
sistent Parameter Values 


The analyses, discussed in previous sections, dealt with determining whether there are 
differences in the parameters of the life equation between companies. This section investigates 
(a) whether three different bearing types made by company B have common values for all the 
parameters in the life equation, and (b) whether the exponents p calculated for each bearing 
tvpe are consistent with the value of p=3. 

The analysis for (a) is similar to the analysis made in section 6.3 of this appendix; i. e., 
separate life equations were fitted to each bearing type and the resulting residual sum of 
squares was compared with the residual sum of squares arising from fitting a single equation 
to all data from company B, irrespective of bearing type. The variance ratio for statistically 
testing (a) is given by eq (C27). 

The analysis for (b) was governed by the following considerations. If the true (or popula- 
tion) value of the exponent p is p=3, regardless of bearing type, then the estimates for p 
obtained by fitting a separate life equation to each bearing type should not differ from p=3 by 
more than the dispersion inherent in the endurance lives of the bearings. The agreement of 
the values of p estimated for each bearing type with p=3 is tested for statistical significance 
by the variance ratio (C28). 

Mathematical formulation. The 148 test groups from company B can be divided into 
three bearing types corresponding to 37 groups for B-1, 94 groups for B—2, and 17 groups for 
B-—3-type bearings. Let these types be denoted by v=1, 2, 3, respectively. Also define 





. 
" 
Af? => @. 102% 1, 7=0, 1, 2,3; v=1, 2, 3 
a=! q 
(C24) 
ne 
Gj? => |) 0,6 2f0 ye =0, 1, 2, 3; vo=1, 2, 3, 
a=! , 


where n, is the number of test groups for bearing type v. Then the normal equations, which 
determine the estimates for the parameters in the logarithmic life equation, are 


, A : i‘ 
> Aj7 oi? =a)” i=0, 1, 2,3; v=1, 2, 3. (C25) 


Thus for each bearing type, the values for the parameters a,(i=0,1,2) and p in the life equation 


can be estimated from the relations 


A 
vr) 
min) _ 95 
a; = ’ 
(pr) 
bs 
a 
A 
p\" — 6s” 
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Therefore the residual sum of squares for the vth bearing type is 


ne 3 
S°= >) Wraia— 201? Gi” v=1, 2, 3, (C26) 
! 


a= i= 


having (n,—4) degrees of freedom. Then the hypothesis of common parameters for the three 
bearing types can be tested by the variance ratio, 
p_(—SY-S®—S)/8 _ 
(SY+ 82+ §®)/136 ’ —" 
having a critical value of F'»;(8, 136) =2.01. 

Since the analysis given in section 6.4 of this appendix reached the conclusion that all 
companies have a common value for p, and since this value (cf. table 2 of main text) with its 
associated uncertainty includes the value p=3 given in [15], it seems also desirable to test a 
second hypothesis that the value of p for each bearing type is consistent with p=3. The 
F-ratio for this hypothesis is given by 


= on as (C28) 
(SY+§%+$®)/136 
where the p‘(v=1,2,3) refer to the estimates of p obtained for each bearing type, and (37 
(yv=1,2.3) is the element occurring in the last row, last column of the inverse matrix to A 7 
(y=1,2.3). The critical value for the variance ratio (C28) is FF», (3, 136) 2.67 

The estimates of a{? (¢=0, 1, 2) and p™ obtained from the solutions of the normal 
equations (C25) are summarized in tables 7 and 4, respectively, in the main text. The caleu- 
lations for the variance ratio (C27) are summarized in table C-2. 

Corresponding to the hypothesis that the three bearing types have the same parameters 
in the life equation, the calculated variance ratios (C27) yield 


a é 
C29) 


TABLE C 2. Synopsis of calculations to determine whether all hear ing type s have common parameters in life 
quation (company B only 


L / 
Degrees 
Type of bearing of 
freedom So Mean S Mean 
square square 
All types combined *144 *285. 121962 F201. 935797 
B-1 33 68. 558208 33. 016973 
B-2 0 174. 325939 131. 142017 
B-3 13 11. 537480 8. 505638 
Sum 136 )«—-254. 421716 1. S707 172. 664628 1. 2006 
Difference Ss 30. 700246 3. 8375 29, 271169 3. 6580 
, - , &ea6a = _ 3.6589 
From eq (C27 fk = z 2.05 f : .=2.88 
l S/O, ] 2606 


*From table C-1 
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Because the critical F value is F'o; (8, 136)=2.01, both the Ly and Lg calculated variance 
ratios are statistically significant. Therefore, one could conclude from the above F-ratios that 
the data support the hypothesis that parameters do differ between bearing types. However, this 
does not exclude the possibility that the values for the exponent p are consistent with the value 
p=—3. Substituting the appropriate quantities in eq (C28) results in 


(C30) 


where 
Cy? =0.53294 


C3 =0.107395 


C2 =0.396744. 


5 


Thus, as both of the above calculated F-ratios are less than the critical value, F 9; (3,136) 

2.67, the conclusion can be made that the values of p are consistent with p=38 for different 

bearing types, although possibly some (or all) of the other parameters in the life formula (i. e., 
a,(?=0,1,2) may differ among the different bearing types. 

The values for a)? ({=0,1,2; e=1,2,3) arising from the analysis by ball-bearing types 
have very large confidence limits (intervals of uncertainty). This is mainly due to the fact that 
an analysis restricted to one bearing type is essentially an analysis on bearings having almost 
the same values for Zand D,. In order to estimate the aj with good precision, it is necessary 
to have results for bearings having wide variations with respect to Z and D,. Thus the esti- 
mates for a; based on all bearing types for company B (table 5 or table 6) have substantially 
smaller confidence intervals as compared to the intervals based only on a single beariag type. 


6.6. Determination and Analysis of f, Based on the Parameter Values 
a, = 2/3, a,=1.8, and p=3 


The values for the parameters a,, ad, and p given in [15] are a,—2/3, aa=1.8, and p=3. 
If these parameter values are valid for the data at hand, then more precise estimates for the 
“workmanship” parameter do (or f,) can be made for each company or bearing type. These 
will generally have better precision compared to the estimates of ad made when the other 
parameters in the life equation are simultaneously estimated along with a. This section 
considers the problem of verifying whether the parameter values, given above, are valid for the 
given data, and for those cases where this is true, estimates of do (or f,.) are obtained assuming 
these values for the other parameters. 

The procedure for determining whether the values a,;=2/3,a,=1.8, and p=3 are valid for a 
given classification of the data (with respect to a company or bearing type) is to fit the data to 
the life equation, using the assumed values for a,, a2, and p. Thus there is only one unknown 
parameter, do, in the life equation to be estimated. Then the resultant residual sum of squares, 
denoted by RP (eq (C37)), can be calculated having (n—1) degrees of freedom. Alternatively, 
the life equation can be fitted to the data such that a// the unknown parameters are simultane- 
ously estimated. The residual sum of squares from this latter fit, S, will have (n—4) degrees 
of freedom. Then, if the above parameter values are not consistent with the given data, R 
will be appreciably larger than S. The variance ratio (C38) is used to determine if the differ- 
ence between these two residual sums of squares, i. e., (2—S) having |(n—1)—(n—4)]=3 
degrees of freedom, is statistically significant. 

Mathematical formulation. Let n be the number of test groups within a particular classifi- 
cation (either by company or bearing type). Then, assuming the values a,;=2/3, a.=1.8, and 
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p=3, the logarithmic life equation for the ath test group can be written as 


Y.=3{dot+7Ze a=1,2,...,%, (C31) 


») 
La E Pia T 1.8 L2a7 ram } (¢ 132) 
» 


The resulting normal equation for estimating dp is 


where 


w,? = 
>> eee... >" w 


hed ata 
A a=l 0” a=! Vee 
ado , (¢ 33) 
n 
—_ 
> j Wa 
a 1 


which also can be written as a function of the sums of cross products g;, d;), &. g., 


> 


A l Yo 2 ~ 
ado a { _ 3 (lo, 7 L.S doe (lo )\. (¢ 34) 


. . » A . 
Hence the estimated variance of do Is 


N 


. A — 
Variance (dp) ’ (C35) 
Vo 
where 
n y e A 
91 S\ow (: ~—J ) — (yz 
_. & b 3 (C36) 
i ] Mi l 


The residual sum of squares ?, having (n—1) degrees of freedom, can also be written as a fune- 


tion of the sums of cross products, 


Ew a 2 } 4 -_ 
R o{5 p WaYa—o $(y—1.2 gots Ist Q @),;+-3.24 dso+a 2.4 dias ay 3.06 dos toot . 
. a 1 ‘ ) ‘ ») 
(C37) 


In the analyses made in the preceding sections, the ball-bearing data have been analyzed 
with respect to individual companies or bearing types. It thus seems desirable to determine 
whether the data within these classifications support the hypothesis that @,—2/5, a,—1.8, and 


p=3. The variance ratio used to test this hypothesis ts 


(R S) > 


F (C38 
S/(n—4) ' 
where, for testing within companies, 
S, for company A ) 
S S, for company B t defined by eq (C13) and given in table C-1, 


S; for company C } 
and for testing within bearing types, a{? is replaced by Af”, 
S” for B-1 
S S® for B-2 defined by eq (C26) and given in table C-2, 
S® for B-3 
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and BR (eq (C37)) refers to the calculated residual sum of squares within the particular 
classification. 

Table C—3 summarizes with respect to companies A, B, and C the values of @, R, the 
calculated F-ratio (eq C38), and the critical F value. Table C-4 summarizes the same quan- 
tities for the B-1, B—2, and B—3-type bearings made by company B. 

The results of this analysis show that the values a;—2/3, a,=1.8, and p=3 are consistent 
for rating life Ly) data with respect to each of the three companies. However, for median 
life Leo, these assumed parameter values are consistent only for the data from company A. 

The results for the same analysis made on the three different bearing types, indicate 
that the assumed parameter values are consistent only for the B-1 and B—2-type bearings 
with respect to median life Lso. 

It is interesting to note that the analysis for company B (ignoring bearing types) showed 
that the assumed values for the parameters are consistent with the Ly) data. However, a 
finer analysis by bearing type revealed that these values are not valid for the B—3-type bear- 
This apparent inconsistency stems from the fact that the analysis for company B, 
taken as a whole, is dominated by those bearing types having the larger number of test 
i. e., B-1 and B-2, and for these types the parameter values were found to be 


ings. 


groups, 
consistent with the Ly) data. 

The estimates for /,, assuming a@,;=2/3, d2=1.8, and p=3, are summarized in table 8 of 
the main text for rating life Ly. This summary also includes the value of f, for company 
D computed from only three test groups. (Because of the small number of test groups, it 
was not possible to verify whether the assumed parameter values are valid for these data.) 


Tasie C-3. Summary of computations (by companies) to test hypothesis that data are consistent with assumed 


values a,;= 2/3, a,a=1.8, p=3 
Lio Lo 
Company Critical 
PF 
4)=I1n f, R F 4)=In f, R F 
\ 8. 4205 102. 539286 *—0. 06 8. 9382 94. 367916 0. 79 2. 80 
B 8. 5021 287. 613198 0. 42 8. 9254 217. 560672 3. 71 2. 65 
5 8. 1001 8. 375562 0. 96 8. 5832 8, 488584 9, 92 3. 59 


*Negative value due to rounding in calculations 


Taste C-4. Summary of computations (by bearing type) for company B to test hypothesis that data are consistent 


with assumed values a;= 2/3, ag=1.8, p=8 
Lio Liso 
Type __—_| Critical 
diy In f, R F dy=In f, R F 
B-1 &. 4575 73. 821978 0. S4 8. 8636 $4. 883126 3. 95 2. 86 
B-2 8. 5236 I81. 128312 1.17 8. 9482 141. 731604 2. 42 2. 70 
B-3 8. 5203 27. 014247 5. 81 9. 0022 18. 905517 5. 21 3. 24 
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